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The Mechanics of Moving Vertical 
Fluidized Systems’ 


III. Application to Cocurrent Countergravity Flow 


D. L. STRUVE?, L. LAPIDUS* and J. C. ELGIN: 


Previous studies in this laboratory have verified 
the assumption of a single fundamental relation be- 
tween void fraction and slip velocity for counter- 
current flow and batch fluidization of solid rigid 
particles in water. The present investigation was 
undertaken to test the applicability of a generalized 
theoretical analysis to the case of cocurrent counter- 
gravity flow. 

The cocurrent countergravity flow and batch 
fluidization of two ideal fluidized systems were 
studied experimentally in a 1 in. diameter column. 
The systems were 0.0184 in. diameter glass spheres- 
water and 0.00396 in. diameter glass spheres-water. 
Holdup, pressure drop, and fluid and particle velo- 
cities were measured for the cocurrent counter- 
gravity flow experiments. Holdup and fluid velocity 
were measured for the batch fluidization runs. 


Calculated slip velocities for cocurrent counter- 
gravity flow and batch fluidization gave a single 
curve when plotted against holdup. It is therefore 
concluded that the experimental results support the 
validity of the basic assumption for the cocurrent 
countergravity flow of liquid-solid systems. 


The experimental data for cocurrent counter- 
gravity flow are in excellent agreement with the 
operating diagram (holdup as a function of fluid 
and particle velocities) for cocurrent countergravity 
flow determined from the holdup-slip velocity re- 
lationship obtained from the batch fluidization ex- 
periments. Stable operation in the cocurrent coun- 
ter gravity flow region where the fluid velocity is 
less than the single particle transport velocity is 
possible provided there is mechanical support for 
the solids and the solids flow rate is controlled at 
the bottom of the column. The application of the 
theory to the prediction of behaviour in compound 
cocurrent countergravity flow fluidized systems is 
demonstrated. 


1Manuscript received May 2, 1958. 

2Esso Research and Engineering Company, Linden, N.]. . 
3Department of Chemical Engineering, Princeton University, Princeton, 
I 


4De partme nt of Chemical Engineering, Princeton University, Princeton, 
N.J. 


This article is based on a paper which was presented at the Joint 
ae h.E.-C.1.C. Chemical Engineering Conference, Montreal, Que., April 
3, 1958. 
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ry two prior publications (1,*) this laboratory has present- 
ed a theoretical analysis of vertical moving fluidized 
systems. The theory considered the ideal or model case of 
a system of rigid, ‘solid spheres suspended by a liquid in 
a vertical tower. It assumes as a basic premise that a 
single fundamental relation exists between the holdup (or 
void fraction) and the slip velocity between particle and 
fluid irrespective of the relative direction of flow of fluid 
and particles relative to each other and relative to the 
walls. This holdup- slip velocity relationship is assumed 
to be unique for a given fluid- -particle system and there- 
fore to be the same for any mode of operation. The 
reader is referred to the first reference for details of this 
development), The second paper(?) illustrated the ap- 
plication of the theory to the case of free countercurrent 
operation of glass spheres or lead shot fluidized by water. 
The results showed excellent agreement between theory 
and experimental data. 


In the present paper the theory is applied to the re- 
sults of an experimental study of free and restrained 
cocurrent countergravity flow of glass spheres and water. 
Figure 1 presents a schematic diagram which illustrates 
the operational arrangement for these types of fluidized 
systems and also the portion of the generalized operating 
diagram which is of interest. This latter diagram plots 
the holdup (1—e) vs. the superficial velocity V; with 
superficial solids velocity Vj,’ as parameter), For co- 
current countergravity flow, with particles denser than 
the fluid, both the fluid ‘ead particle velocities are 
directed upward and therefore under the convention 
adopted are positive in sign. When there is no mechani- 
cal support beneath the particles, fluid velocities greater 
than the single particle transport velocity are required 
for cocurrent countergravity flow. This type of opera- 
tion is called free cocurrent countergravity ‘flow or trans- 
port. A mechanical support or restraint beneath the 
particles is a necessary condition for cocurrent counter- 
gravity flow at fluid velocities less than the single particle 
transport velocity. This type of operation is called bot- 
tom restrained cocurrent countergravity flow. 

The cocurrent countergravity flow and batch fluid- 
ization of two ideal fluidized systems were studied ex- 
perimentally in a 1 inch diameter column. The systems 
were 0.0184 in. diameter glass spheres-water and 0.00396 
in. diameter glass spheres-water. Holdup, pressure drop 
and fluid and particle velocities were measured for the 


141 








| 

| ° 
o 2.6 

| 50 

! aes 
° ° 

° Co 

| vet @Xtrecyes® 

, ———— Screen 

| 

| 

| 

| 

| 

| 

| 







Transport Velocity of Single Particle 


| 
| 
| 
| 
| 
| 
| 
i 
| 
| 
| 
| 
| 
| 
| 
| 
| 







39,0 
0,0 
ol, 
Vy = + 2 — ° 
y= + 
ssn ee sig aaa es aia es a 





FREE COCURRENT _ 
COUNTERGRAVITY FLOW ~ 
(Transport) 





FLUID VELOCITY, Vy 


(/-€) 3 X 
\ 

\ mm \ 
> BOTTOM RESTRAINED 
I COCURRENT COUNTER-~ 
~ GRAVITY FLOW 
a, + 
4 % 

2 ‘, 
0 
Figure 1 — Schematic diagram of various fluidized 


cocurrent countergravity flow experiments. Holdup and 
fluid velocity were measured for the batch fluidization 
runs. The experimental results support the validity of the 
basic assumption of a unique holdup- slip velocity re- 
lationship. In addition the data are in excellent agreement 
with the operating diagram predicted on the basis of only 
the batch fluidization ex ‘periments. The application of 
the theory to the prediction of behavior in compound 
cocurrent countergravity flow fluidized systems is also 
demonstrated 


Previous studies of transport 


Pneumatic transport of solids has been studied by 
several investigators in order to measure and correlate 
pressure drops in vertical and horizontal systems. The in- 
vestigations of Gasterstadt(*) and Cramp and Priestly“) 
on grain-air systems, Belden and Kassel “), and Farbar (*) 
on catalyst-air systems, Hariu and Molstad) on sand 
and cataly st-air systems, and Vogt and White ‘*) on sand, 
steel-shot, clover seed, and wheat-air systems report pres- 
sure drop measurements for transport in vertical columns. 
Correlations of the total pressure drop were developed 
which treated the total pressure drop as the sum of in- 
dividual pressure drops to pipe wall-fluid friction, pipe 
wall-solid friction, acceleration of the fluid, acceleration 
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systems and operational diagram for the case of transport. 


of the solid, and the static pressure drop or weight of 
solids and fluid in the column. 

In all reported cases, appreciable pipe wall-fluid fric- 
tion and acceleration of the solids caused the total pres- 
sure drop to be much higher than the static pressure 
drop. The static pressure drops given by Belden and 
Kassel may be in error as the authors assumed that for 
a given fluidized system the slip velocity was constant 
and independent of holdup. Cramp and Priestly neglected 
the increase in fluid velocity above the superficial fluid 
velocity due to the decreased area available for flow in 
the presence of particles. In all cases where the holdup 
was given, or could be calculated from the reported data, 
there was no simple relationship between slip velocity 
and holdup as the solids had not been allowed to ac- 
celerate to their final steady-state velocity. 

Zenz‘) pointed out the relationship between total 
pressure drops for batch fluidization and those for verti- 
cal transport in air-solid systems. The available data did 
not enable him to investigate the relationship between 
static pressure drop or holdup for batch fluidization and 
vertical transport. 

Bhattacharya znd Roy “') and Mertes and Rhodes ‘") 
present data for the transport of solids by liquids in 
vertical columns. Bhattacharya and Roy transported 
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Fischer-Tropsch catalyst with kerosene and water. Total 
pressure drops were measured and a correlation was de- 
veloped following the method of Vogt and White for 
pneumatic transport. They calculated static pressure 
drops using values of holdup equalling the ratio of solids 
volumetric flow rate to total volumetric flow rate. In 
light of the present work the above calculation was er- 
roneous as the holdup is determined by the slip velocity 
and is not equal to the ratio of volumetric flow rates. 
Mertes and Rhodes reported a few data for the vertical 
transport of glass spheres by water. Pressure drops were 
not measured and the range of holdups (0.0726 to 0.184) 
is not sufficiently broad to enable a critical test of the 
theory presented in this paper. 


EXPERIMENTAL 


Materials 


Glass spheres. The glass spheres selected for the 
batch fluidization and cocurrent countergravity flow ex- 
periments were 3M “Superbrite” Glass Beads, Grade Nos. 
107 and 113, manufactured by the Minnesota Mining and 
Manufacturing —— Batch fluidization experiments 
were also done using glass spheres, hereafter referred to 
as No. II spheres, on hand in this laboratory and of un- 
known source. 

The density of each grade of spheres was determined 
by water displacement. Microscopic measurements, with 
a micrometer eyepiece, of the diameter of 100 spheres 
of each grade were made. The diameter of a grade of 
spheres was taken as the arithmetic average of the 100 
diameters for that grade. The No. 107 spheres passed 
U.S. Screen Size No. 35 and were retained on No. 140. 
The No. 113 spheres passed No. 140 and were retained 
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on No. 170 screen. City water was used in all experi- 
ments. The average temperature of the water was 22.8°C. 
The results of the measurements are summarized in 





Table 1. 
TABLE 1 
PHYSICAL PROPERTIES OF GLASS SPHERES 
) : 
Glass | Diameter — 
Sphere | Dp 
No. (in.) | Alb./ oa ait ) 
Mee | 0.0370 159.6 

Mois cveeid 0.0184 155.0 
Bag arse aire | 0.00396 153.1 


Apparatus for cocurrent countergravity 
flow experiments 

A detailed schematic diagram of the cocurrent coun- 
tergravity flow apparatus is shown in Figure 2 and should 
be referred to while following the discussion below. 


Water supply. Water from constant head tank 1 
(25 gal. capacity) flows to a base of the column in suc- 
cession through a Y-strainer, pump 1, a globe valve, 
rotameter 1, 2, or 3, and a gate valve. All pipe is gal- 
vanized iron with brass valves. Make-up water is added 
to constant head tank 1. 


Solids feed to the column. The solids feed line 
from the bottom of hopper 1 to the base of the column 
is 1 in. Pyrex Glass Pipe, for visibility, with the exception 
of short lengths of 1 in. galv anized iron pipe on either 
side of the orifice. Interchangeable plates with circular 
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Figure 2 — Detailed diagram of cocurrent countergravity flow apparatus. 
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orifices from 45/64 to 1 in. diameter are used to control 
the rate of solids flow. The plates are of 4 in brass and 
are inserted between 1 in. flanges in the vertical section 
of the solids feed line, 6 in. above the 45° elbow. The up- 
stream side of the orifice is enlarged to 1 in. diameter. 
The solids enter the base of the column at an angle of 
45° above the horizontal. 

The top of hopper | is sealed from the atmosphere. 
Pressure is applied above the solids by opening the needle 
valve in the by-pass line, which runs from a point fol- 
lowing pump 1 to the top of hopper 1. A vent at the top 
of hopper 1 is provided for the removal of air. Pressure 
above the solids aids in establishing high solids flow rates 
by preventing the backflow of water up the solids feed 
line from the base of the column as well as exerting an 
additional downward force on the solids. 


Column. The base of the column is a special 
tee of Pyrex glass pipe with a sidearm of 1 in. diameter 
at an angle of 45°. A piece of No. 170 screen is inserted 
between the gaskets of the flangles joining the bottom of 
the 2 in. tee to the water supply line. The remainder of 
the column is of 1 in. Pyrex glass pipe and is joined to 
the top of the tee by a special glass pipe reducer with a 
16° taper. From bottom to top the remainder of the 
column is a 3 ft. acceleration section, lower pressure tap, 
lower plug valve, 3 ft. holdup measuring section, upper 
plug valve, upper pressure tap, | ft. straight length, U- 
bend, and glass pipe to the top of hopper 2. The plug 
valves are standard 1 in. Pyrex glass pipe fittings. The 
column is supported by a superstructure of 1} X 1} X 18 
in. angle iron. 

Details of the pressure taps are shown in Figure 3. 
Lines of } in. copper tubing run radially outward from 
each of the three holes to a 4 in. diameter circle of 
copper tubing. A line of n. copper tubing runs from 
the circle to a manometer manifold. The manifold facili- 
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Figure 3 — Details of pressure taps. 


tates the bleeding of air from the manometer lines. The 
bleed lines and lines to the U-tube differential mano- 
meters are also } in. copper tubing. A 6 in. manometer 
filled with mercury and a 127 cm. manometer filled with 
carbon tetrachloride provide good sensitivity over the 
complete range of pressure drops. 

Solids return. The water and solids flow from the 
top of the column to hopper 2, where the water over- 
flows and is returned to constant head tank 1 through 
4 in. leader pipe. The solids settle to the bottom of 
hopper 2. At the end of a run, the solids are conveyed 
to hopper 3 through a Schutte and Koerting Co. 1 in. 
Water Jet Eductor. The solids and their interstital water 
flow from the bottom of hopper 2 to the suction inlet 
of the eductor, where they are mixed with the motive 
water and transported up to hopper 3. 

The motive water for the eductor flows from con- 
stant head tank 2 ag gal. capacity) in succession through 
a Y-strainer, pump 2, and a globe valve to the eductor. 
Pressure gauges 1 and 2 (100 Ib./sq.in. scale) indicate if 
sufficient pressure is being developed to operate the 
eductor. All pipe in this section is galvanized iron. A 
gate valve at the end of the transport line can be closed 
to prevent siphoning back through the eductor to hopper 
2 when transport of the solids to hopper 3 is completed. 
Make-up water is added to constant head tank 2. 


Solids collection and storage. The mixture of 
water and solids from the eductor is collected in hopper 
3 until hopper 2 is empty. The solids settle to the bottom 
of hopper 3 and are drawn off into hopper 1. As the 
solids flow into hopper 1, water is displaced and flows 
through the vent and a 4 in. leader pipe to constant head 
tank 2. The solids settle to the bottom of hopper 1 where 
they are stored, under water, until the next run. When 
the solids have completely settled in hopper 1 the water 
remaining in hopper 3 is returned to constant head tank 
2 via the vent. 

Sketches of hoppers 1, 2 and 3 are shown in Figure 4. 


Procedure for cocurrent countergravity 

flow experiments 
The sequence of events during a run is outlined be- 

low: 

1. The plug valves are opened and pump 1 is turned on. 

2. The needle valve in the by-pass line is partially 
opened. 

3. The desired water flow rate is maintained by the 
globe valve at the rotameter, which measures the 
desired flow rate. 

4. Air is bled from the manometer lines by. means of 
the manifolds and the manometer valves are opened 
so the manometers are operative. 

5. When the manometer has remained steady for 3 or 
4 minutes, the solids-water effluent from the top of 
the column is collected in a graduated cylinder for 
a period of time, measured with a stopwatch. 

6. The manometer is checked to see that it is un- 

changed and the manometer valves are closed. 

The two plug valves are closed simultaneously. 

When the solids have settled, the height of the bed 

formed above the lower plug valve is measured with 

a meter stick. 

8. The volume of solids and water collected in the 
graduated cylinder and the height of the settled bed 
of solids in the cylinder are determined. 

9. The manometer and stopwatch are read. 

10. Pump 2 is turned on and the eductor is operated 
until all the solids in hopper 2 have been transferred 
to hopper 3. The valve at the end of the transport 
line is closed. 


| 
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11. The solids are transferred from hopper 3 to hopper 
1 and the water remaining in hopper 3 is drawn off 
until hopper 3 is only one-third full. The valve at the 
end of the transport line is opened. 


Measurements made during cocurrent 
countergravity flow experiments 

Superficial liquid and particle velocities. Since 
water enters the column with the solids as well as from 
the main water supply, the water flow rate indicated by 
the rotameters is too low. The true water flow rate is 
determined as described below. 

The graduated cylinders are calibrated so the weight 
of solids in the cylinder is known as a function of the 
height of the settled bed of solids measured with a meter 
stick. Knowing the weight of solids collected in unit time 
and the density of the particles, the superficial particle 
velocity V,’ for the 1 in. column can be calculated. The 
total volume of the mixture in the cylinder less the 
volume of solids in the cylinder (weight of solids/density 
of solids) is the volume of water collected in unit time, 
from which the superficial liquid velocity V;’ can be 
calculated. 

Holdup. The section of the column between the 
two plug valves is calibrated so the weight of solids in 
the section is known as a function of the height of the 
settled bed in the column between the plug valves, 
measured with a meter stick. The total volume, deter- 
mined by filling with water, of the holdup measuring 
section is 576 cu. cm. Dividing the weight of solids by 
the density of the solids and the total volume of the hold- 
up measuring section gives the holdup (1—e). 


Apparatus for batch fluidization experiments 


The apparatus used for the cocurrent countergravity 
flow experiments is used for the batch fluidization ex- 
periments with the following modifications: 

1. The solids feed line is disconnected at the base of the 

column and the sidearm of the 2 in. tee is sealed with 

a } in. brass blind flange. 

The plug valves, pressure taps, and 1 ft. length of 

glass pipe are removed and a 3 ft. section of glass 

pipe is added above the two remaining sections of 
glass pipe. 

3. A screen (No. 40 for No. 107 and No. II spheres and 
No. 170 for No. 113 spheres) is inserted between the 
lower and middle 3 ft. sections to act as a support 
for the bed. 

The water from the top of the column flows into 
hopper 2 and overflows to constant head tank 1, as for 
cocurrent countergravity flow. The eductor system is 
not used for batch fluidization. 


m 


Procedure for batch fluidization experiments 
The sequence of events during a run is outlined be- 

low: 

1. A known weight of solids is poured into the column 
from the top. 

2. Pump | is turned on and the water flow rate is in- 
creased until the bed is expanded to the full height 
of the column. 

3. The height of the fluidized bed is measured with a 
meter stick at each of a series of decreasing water 
flow rates. Each flow rate is maintained until the sur- 
face of the bed remains at the same level for 2 or 
3 minutes. 


Measurements made during batch fluidization 
ex periments 


Superficial liquid velocity. The rotameters are cali- 
brated so the true flow rate in gal./min. may be calcu- 
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Figure 4 — Sketch of solids storage hoppers 1, 2 and 3. 


lated from the rotameter scale reading. The true flow 
rate is converted to the superficial liquid velocity V;’ 
for the 1 in. column. 

Holdup. The height of the bed for zero voidage is 
calculated by dividing the weight of solids in the bed by 
the density of the particles and the cross-sectional area 
of the 1 in. column. The holdup (1—e) is given by the 
ratio of this height to the height of the fluidized bed. 


THE HOLDUP-SLIP VELOCITY RELATION FOR 
COCURRENT COUNTERGRAVITY FLOW 


Slip. velocities were calculated for the cocurrent coun- 
tergravity flow experiments by using the measured values 
of holdup and superficial liquid and particle velocities 
in conjunction with Eq. (1) 

re ra 
She gh Fig ac cand nl (1) 
€ l1—e 
Slip velocities were calculated for the batch fluidization 
experiments by using the measured superficial liquid 
velocities, the averages of the holdups measured for dif- 
ferent bed weights at the same superficial liquid velocity 
and Eq. (1). In this latter case, however, V4’ = 0 and 


V, oi Whe alloc yn hen eegee el (2) 


| batch = € 

Holdup is plotted against the calculated slip velocities 
for the cocurrent countergravity flow and batch fluidiza- 
tion of 0.0184 in. glass spheres and water in Figure 5 and 
of 0.00396 in. glass spheres and water in Figure 6. A line 
is drawn through the batch fluidization points. This line 
is extrapolated to the left to the holdup of the packed 
bed (1—e), and is extended to the right to the single 
particle transport velocity at zero holdup. The holdup 
(1—e), of the packed bed formed by slowly decreasing 
the fluid flow rate in a batch fluidization experiment 
and the single particle transport velocity calculated from 
the drag coefficient-Reynolds number correlation ‘'?), are 
given in Table 2 for the three sizes of glass spheres. 

In Figures 5 and 6 the points for cocurrent counter- 
gravity flow fall along the line drawn through the batch 
fluidization points. The agreement between the holdup- 
slip velocity relations for cocurrent countergravity flow 
and for batch fluidization supports the basic assumption 
of the existence of a single holdup-slip velocity relation 
for any vertical fluidized system. This assumption is also 
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Figure 5 — Holdup slip velocity relation for No. 107 glass 
spheres and water from batch fluidization and cocurrent 
countergravity flow experiments. 


TABLE 2 


LIMITs OF THE HoLpup-SLip VELOCITY RELATION 


Single Particle 
Transport Velocity 


Holdup of 
Packed Bed 


Diameter of 
Glass Sphere Dy | 


(in. ) (1 — €)p (ft./min. ) 
0.00396...........| 0.526 1.55 
Ps ivavins oad 0.526 13.9 
aaa 0.517 28.2 


supported by the results of previous studies “': *) of the 
countercurrent flow and batch fluidization of liquid-solid 
fluidized systems. Since the Reynolds numbers based on 
superficial fluid velocity for batch fluidization are less 
than those for free cocurrent counter rgravity flow, the 
results of the present study indicate there is no effect of 
Reynolds number based on superficial fluid velocity on 
the holdup-slip velocity relation for a given liquid- -solid 
system. The experimental results also show no effect of 
Rey nolds number based on superficial particle velocity 
on the holdup-slip velocity relation for a given liquid- 
solid system. 

The above conclusions for cocurrent countergravity 
flow and batch fluidization and previous conclusions for 
countercurrent flow and batch fluidization are drawn 
from experimental results obtained with two or more 
fluidized systems of different physical properties. Taken 
together, these conclusions support the generalization 
that for any vertical fluidized system a single fundamental 
relation exists between the holdup (or void fraction) 
and the slip velocity. 


Statistical analysis of the experimental results 


Statistical analysis of the experimental results was 
pursued in order to obtain further support for the exist- 
ence of a single relation between holdup and slip velocity 
for the cocurrent countergravity flow and batch fluidiza- 
tion of fluidized systems. 


I — Batch fluidization results. The results of the 
batch fluidization experiments indicate that for constant 
liquid flow rates the holdup varies with the weight of 
bed. To obtain a single holdup-slip velocity curve from 
the batch fluidization results for a system, the holdups 
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Figure 6 — Holdup slip velocity relation for No. 113 


glass spheres and water from batch fluidization and co- 
current countergravity flow experiments. 


at a constant liquid flow rate were averaged and the 
average holdup used to calculate the slip velocity from 
Ba. (2). 

The spread of the actual batch fluidization data about 
the single smooth curve representing them (as plotted 
in Figures 5 and 6) is indicated in Table 3. 

For each system the constants of an equation of the 
form: 

V; =a+bi(1 — €) + be(1 — €)?..... ee 
were determined by curvilinear regression of the raw 
batch fluidization data. The values of the constants and 
their standard deviations are given in Table 4. It is 
demonstrated in Figures 7 and 8, where points calculated 
from Eq. (3) are compared wtih the single smooth curve 
representing the average experimental batch fluidization 
data for the 0.0184 in. and 0.00396 in. spheres, that an 
empirical equation of this form fits the batch fluidization 
data. 


Il — Cocurrent countergravity flow results. The 
scatter of the points for cocurrent countergravity flow 
on Figures 5 and 6 is the result of two factors: experi- 
mental error inherent in the method used to measure 
holdup and fluid and particle velocities, and failure to 
attain steady-state behavior. 

The experimental error of the results is estimated to 
be: 

(1—e) accurate to the third decimal place 

Vi + 1.5% 

Va = 1% 

As there is no significant error in the holdup (or void 
fraction), the scatter of the points on the holdup-slip 
velocity plot arises from the error in the calculated slip 
velocity. The uncertainty in the slip velocity is due to the 
errors in the fluid and particle velocities which are used 
in the calculation of the slip velocity from Eq. (1). Sub- 
stitution of the velocities and their errors in Eq. (1) 
gives 
V (Ve) 1.5%) (Va + 1%) 

. € (i-¢«) 
(Ve + 1.5%) — (Va + 1%)...... 
and the error in the slip velocity depends on the magni- 
tude of the fluid and particle velocities. In this equation 
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V, = average fluid velocity, ft./min. and V, = average 
solids velocity, ft./min. The results of the esteslation of 
the uncertainty in the slip velocity due to experimental 
error for four of the more w idely scattered points from 
Figure 6 are given in Table 5. The values of slip velocity 
from batch fluidization experiments at the same holdup 
are given for comparison. 


Fluidized systems require time to attain their steady- 
state conditions. Operating experience gained in this 
work showed that often 10 to 15 minutes are required 
until stable operation occurs. The experimental criteria 
for the attainment of steady-state is the constancy of the 
manometer reading for at least two or three minutes. 

The time available for a run was limited by the capa- 
city of the solids feed hopper. Greater solids weight 
flow rates were required for the No. 107 spheres than for 
the No. 113 spheres. As the solids feed hopper had the 
same weight capacity for both size spheres, more time 
was available for a run with the No. 107 spheres than 


TABLE 3 


SPREAD OF BATCH FLUIDIZATION DATA 








0.0184 in. Spheres 0.00396 in. Spheres 
(1 —¢) Vv. |a-e \. 
(ft./min. ) (ft./min. ) 
High bed weight | 0.510 | 1.05 | 0.488 0.0938 
Low bed weight 0.525 | 1.08 0.476 0.0916 
Average for line 0.517 | 1.07 0.481 0.0925 
High bed weight | 0.349 | 3.09 | 0.340 | 0.306 
Low bed weight 0.376 | 3.22 0.330 0.301 
Average for line 0.363 | 3.16 | 0.335 0.304 
High bed weight | 0.246 5.09 | 0.267 0.491 
Low bed weight | 0.264 Sues | 0.256 | 0.484 
Average for line 0.253 5.14 0.261 | 0.487 
High bed weight 0.070 | 9.06 0.101 0.973 
Low bed weight 0.070 9.06 0.098 0.970 
0.100 0.972 


Average for line 0.070 9.06 
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Figure 7 — Comparison of empirical batch fluidization 


equation fitted by curvilinear regression with experimental 
results for batch fluidization of No. 107 glass spheres 
and water. 
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with the No. 113 spheres. In some runs with the No. 107 
spheres steady-state conditions had not quite been 
reached when the solids supply ran low, and the data for 
the run were taken with the manometer reading slowly 
changing. If the manometer reading was changing rapidly 
the run was discarded. 

The greater spread for the No. 107 spheres (Figure 5) 
than for the No. 113 spheres (Figure 6) is a reflection of 
the departure from steady-state behavior associated with 
the No. 107 spheres. 

For each system, the constants of Eq. (3) were also 
determined by curvilinear regression of the cocurrent 
countergravity flow data. The values of the constants 
and their standard deviations are given in Table 4. 


Comparison of batch fluidization and cocurrent 
countergravity flow results 

Examination of Figures 5 and 6 in light of the spread 
of the batch fluidization data (Table 3) and a degree of 
uncertainty in the calculated slip velocity for cocurrent 
countergravity flow (Table 5) further supports the 
“graphic” demonstration of the existance of a single 
holdup-slip velocity relation for batch fluidization and 
cocurrent countergravity flow. The greater scatter for 
the No. 107 spheres (Figure 5) is attributed to deviations 
from steady-state behavior which are not a factor for 
the No. 113 spheres (Figure 6). 

“Mathematical” confirmation of the existence of a 
single holdup-slip velocity relation for cocurrent coun- 
tergravity flow and batch fluidization may be obtained 
by comparing the values of the constants b, and b, of 
Eq. (3) calculated from the batch fluidization data ‘and 
from the cocurrent countergravity flow data. The 
hypotheses: 


bip = bic 


bow = Dee 


are made for both systems studied. A 10% risk (a = 10%) 
~ acceptable in the rejection of any of the hy potheses. 
Values of the “Student” t calculated from the equation: 


beocur — Dpatch 


“Student” t = t, = ————___——_....... . (5) 


S*b — S*p 


cocur batch 


0.00396 IN. GLASS SPHERES-WATER 


BED (I-€)p 


co 








— AVERAGE EXPERIMENTAL VALUES 
FROM BATCH FLUIDIZATION 


© CALCULATED FROM 
Vg = 1.386 —4.294( I-€ )+3.245(I-«)* 


HOLDUP (I-€) 
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Figure 8 — Comparison of empirical batch fluidization 
equation fitted by curvilinear regression with experimental 
results for batch fluidization of No. 113 glass spheres 
and water. 
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RESULTS OF 


COUNTERGRAVITY FLow Data FoR No. 107 AND 113 GLAss SPHERES AND WATER 
V; =a+bi (1 —€) + be (1 — €) 
No. 107 Spheres No. 113 Spheres 
D, = 0.0184 pa = 155.0 Dp = 0.00396 Pa 153.1 
(in. ) (Ib. /cu.ft.) (in. ) (lb./cu.ft. ) 
Batch Cocurrent Batch Cocurrent 
Fluidization Countergravity Fluidization Countergravity 
Flow Flow 
Bice skank + eae 10.64 9.513 1.386 1.316 
ei * 25.31 —15.70 —4.294 —3.729 
be ..43 2.535 3.245 2:35 
Spi. 0.4868 8.866 0.2190 0.3009 
Spo. 0.7469 21.52 0.3755 0.5317 
Degrees of Freedom...... 160 42 69 27 
COMPARISON OF CONSTANTS FROM BATCH FLUIDIZATION AND COCURRENT COUNTERGRAVITY FLOW 
Deo. ur ~ Dpat: h 
t. = 
[S*beox ur S*bbatch}! 2 
No. 107 Spheres No. 113 Spheres 
b; b b; be 
t, from above equation and data.... 1.08 0.729 1-33 1.16 
tioe, from tables “» ; 1.68 1.68 1.70 1.70 
TABLE 5 
UNCERTAINTY IN SLIP VELOCITY DUE TO EXPERIMENTAL ERROR FOR COCURRENT COUNTERGRAVITY FLOW 
V;’ Va’ Ve _ Va _ Vs Veb 
(1 — €) (ft./min. ) (ft./min.) (ft./min. ) (ft./min. ) (ft./min.) (ft./min. ) 
0.104 | 3.36 0.282 StS 0.056 zua8 0.027 1.04+0.083 0.96 
0.276 2.80 0.972 3.88 0.058 3.33 0.035 0.35+0.093 0.45 
0.411 0.384 0.169 0.65 0.010 0.41 0.004 0.24+0.014 0.16 
0.519 2.07 2.22 1.30 0.065 4.28 0.043 0.02+0.108 0.075 
| 
are given in Table 4. The corresponding values of t,,,. by the results of a statistical analysis of the experimental 


are also listed. In all cases, the value of t calculated from 
Eq. (5) is less than t,, %. Therefore, the hypotheses pre- 
sented above for both systems can not be rejected at the 
10% risk level. A statistical analysis of the experimental 
data thus lends additional support to the assumption of 
the existence of a single holdup-slip velocity relation for 
batch fluidization and cocurrent countergravity flow. 
The holdup-slip velocity relations calculated from 
Eq. (3) for batch fluidization and cocurrent counter- 
gravity flow are compared in Figure 9 for the No. 107 
spheres and in Figure 10 for the No. 113 spheres. The 
agreement between the curves for batch fluidization and 
cocurrent countergravity flow is excellent in both cases. 
In summary, the conclusions, stated after cursory 
examination of the results presented in Figures 5 and 6, 
supporting the assumption of the existence of a single 
fundamental relation between holdup and slip velocity 
for any vertical fluidized system are further strengthened 
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TABLE 4 
CURVILINEAR REGRESSION AND STATISTICAL 


ANALYSIS OF 


THE BATCH FLUIDIZATION AND COCURRENT 


data. 


HOLDUP AS A FUNCTION OF FLUID AND PARTICLE 
VELOCITIES IN COCURRENT COUNTERGRAVITY 
FLOW 


The theoretical curves of holdup as a function of fluid 
velocity for constant values of particle velocity for co- 
current countergravity flow are shown in Figure 11 for 
the No. 107 glass spheres and water in Figure 12 for the 
No. 113 glass spheres and water. 

The characteristics of the cocurrent countergravity 
flow apparatus used in this work did not permit opera- 
tion at a series of liquid velocities for a constant particle 
velocity. However, by coincidence, certain particle velo- 
cities were repeated and these experimental points are 
shown in Figures 11 and 12. A few isolated experimental 
points are also shown to indicate the range of variables 
covered in the experimental work. 
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Figure 9 — Holdup slip velocity relations for No. 107 


glass spheres and water from curvilinear regression of 
batch fluidization and cocurrent countergravity flow 
experimental results. 
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Figure 11 — Cocurrent countergravity flow region of the 


generalized operating diagram for No. 107 glass spheres 
and water. 


The agreement between theory and experiment is ex- 
cellent for both sizes of spheres. The lines through the 
experimental points for cocurrent countergrav ity flow 
have the proper position and curvature with respect to 
the theoretical lines predicted from the batch fluidization 
data. The theory, as applied to cocurrent countergravity 
flow, demands that for decreasing liquid velocity for any 
constant particle velocity there is a smooth transition 
from zero holdup to the holdup of a packed bed (1 —e),. 
The lines through the experimental points for cocurrent 
countergravity flow in Figures 11 and 12 support the 
above statement. The two ideal fluidized systems studied 
in this work exhibited no discontinuities at intermediate 
values of holdup due to slugging and choking as de- 
scribed by Zenz ) for non-ideal gas-solid transport. 


FREE AND BOTTOM RESTRAINED COCURRENT 
COUNTERGRAVITY FLOW 
For vertical fluidized systems the theory predicts that 
stable cocurrent countergravity flow is possible with 
superficial fluid velocities less ‘than, as well as greater 
than, the single particle transport velocity. However, in 
order to achieve stable cocurrent countergravity flow 
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Figure 10 — Holdup slip velocity relations for No. 113 


glass spheres and water from curvilinear regression of 
batch fluidization and cocurrent countergravity flow 
experimental results. 
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Figure 12 — Cocurrent countergravity flow region of the 
generalized operating diagram for No. 113 glass spheres 
and water. 
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with fluid velocities below the single particle transport 
velocity it is essential that a mechanical support or re- 
straint for the particles be provided at the bottom of the 
column. It is also essential that the particles be removed 
from the top of the column at the same rate at which 
they are fed to the bottom. The necessity of mechanical 
support for the particles at fluid velocities below the 
single particle transport velocity provides the basis for 
distinguishing the following two types of cocurrent 
countergravity flow: free (no mechanical support re- 
quired) and bottom restrained (mechanical support re- 
quired beneath the particles). Free cocurrent counter- 
gravity flow is also called transport. 

The presence or absence of a mechanical support at 
the base of a column should not influence the stability 
of a system in free cocurrent countergravity flow. The 
theory predicts, as shown by the theoretical lines on 
Figures 11 and 12, that, provided a mechanical support 
is present at the bottom of the column, there is a smooth 
transition between bottom restrained and free cocurrent 
countergravity flow. The apparatus used in this work 
was designed to permit both bottom restrained and free 
cocurrent countergravity flow operation. The experi- 
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mental points on Figures 1] and 12 confirm the predic- 

tion of a smooth transition between free and bottom 

restrained cocurrent countergravity flow when the equip- 
ment is properly designed with a restrained support for 
the particles. 

The requirement of mechanical support or restraint 
for the particles at the bottom of a column when fluid 
velocities in cocurrent countergravity flow are below 
the single particle transport velocity is the basis for de- 
signating such systems as bottom restrained. The me- 
chanical support may be either a screen or a more dense 
fluidized system. 

As a second requirement for stable operation of bot- 
tom restrained cocurrent countergravity flow systems 
provision must be made for the particles to leave the top 
of the column at a rate equal to that at which they are 
fed into the bottom terminal. If the particles can not 
flow from the top terminal as rapidly as they enter the 
the bottom, a fully packed bed will form at the top. This 
bed will grow toward the feed entrance until the entire 
column becomes choked with solids and operation ceases. 

By starting with a batch fluidized bed (V4 = 0) and 
then commencing to feed solids to the Sabet: the bed 
is caused to expand and overflow the top of the column. 
A bottom restrained cocurrent countergravity flow fluid- 
ized system is thus obtained. This system has a higher 
holdup than that in a batch fluidized bed at the same 
superficial liquid velocity. 

The apparatus designed for this work actually op- 
erated as a compound column. The upper 1 in. section 
was joined to the top of a short 2 in. section by a tapered 
section. Solids were fed to the 2 in. section, beneath 
which there was a screen. The superficial liquid velocity 
in the 2 in. section was one-fourth that of the superficial 
liquid velocity V;,’ in the | in. section, With this design, 
therefore, three types of compound fluidized systems 
were possible: 

1. Free transport in both the | in. and 2 in. sections. 

2. Free transport in the | in. section and bottom re- 
strained cocurrent countergravity flow in the 2 in. 
section, 

3. Bottom restrained cocurrent countergravity 
both the | in. and 2 in. sections. 

For free transport in both sections the superficial 
liquid velocity V, in the | in, section was greater than 
four times the single particle transport velocity. No 
mechanical support was then needed. 

For free transport in the | in. section and bottom 
restrained cocurrent countergravity flow in the 2 in. 
section the superficial liquid velocity V,’ in the 1 in. sec- 
tion must exceed the single particle transport velocity 
and be less than four times the single transport velocity. 
The superficial particle velocity V,’ in the 1 in. section 
must not be greater than four times the critical value of 
the superficial particle velocity (the value of V,’ for 
V; = single particle transport velocity and (1—e),). 
Mechanical support is needed for the 2 in. section. The 
screen served this purpose. 

For bottom restrained cocurrent countergravity flow 
in both sections the superficial liquid velocity V;’ in the 
1 in. section must be at least four times the fluidization 
velocity for a particle velocity one-fourth of the super- 
ficial particle velocity V,’ in the 1 in. section. The super- 
ficial particle velocity in the 1 * section cannot exceed 
the critical value. Mechanical support is needed for both 
sections. The screen serves for the 2 in. section. How- 
ever, a screen cannot be placed between the 2 in. and 
1 in. sections without I & g the solids flow. Fortun- 
ately, the more dense (higher holdup) fluidized system 
in the 2 in, section provides the necessary mechanical 


flow in 
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support for the less dense fluidized system in the 1 in. 
section. 

The experimental data, taken in the | in. section, show 
that the compound column used in this work gave stable 
operation for any of the three combinations of fluidized 
systems outlined above. The experimental results for 
three runs corresponding to operation in each of the 
three combinations of fluidized systems are shown in 
Table 6. 

TABLE 6 
COMBINATIONS OF COCURRENT COUNTERGRAVITY FLOW 
FLUIDIZED SYSTEMS IN A 1 IN. ABOVE 2 IN. COMPOUND VERTICAL 
Cotumn. No. 113 GLass SPHERES AND WATER 


Single Particle Transport Velocity = 1.55 ft./min. 
V;’ V;’ | Va' (1 — e€) 
in the in the in the in the 
1 in. 2 in. 1 in. 1 in 


section | section section section 


(ft./min.)/(ft./ min. )/(ft./min. ) 





Free transport in | 

both sections 6.98 | 1.74 | 0.994 0.137 

Free transport above | 

bottom restrained 

cocurrent counter- — | 

gravity flow 2.80 0.700 0.972 | 0.276 

Bottom restrained 

cocurrent counter- 

gravity flow in 

both sections 1.0 2 «| 6.275 0.480 0.343 
| 





The experimental data obtained in this work demon- 
strate, for the first time, the stable operation of bottom 
restrained cocurrent countergravity flow fluidized sys- 
tems as predicted by the theory. ; 


PRESSURE DROP IN COCURRENT 
COUNTERGRAVITY FLOW 


Previous workers have assumed the total pressure 
drop AP to be the sum of the static pressure drop AP,, 
fluid-wall friction AP,, solids-wall friction, a solids ac- 
celeration loss, and a fluid acceleration loss. 

Total pressure drops AP over a 5.42 ft. length of 1 in. 
column were measured by means of suitable manometers 
for the cocurrent countergravity flow runs. The esti- 
mated error for the total pressure drop (due to reading 
of the manometers) increased from 0 to £1% as AP 
increased from 0 to 250 Ib./sq.ft. The fluid-wall friction 
pressure drop AP; can be estimated from the standard 
friction factor correlation for smooth tubes using the 
actual fluid velocity (V;= V;,’/e) as it is representative 
of the average fluid velocity in the fluidized system. The 
corrected pressure drop AP. can be calculated from: 

Al. = AP — AP; 
and the static pressure drop calculated from the Eq. 
AP, = Z(1—€) (pa—pe)s using the experimental values 
of holdup. 

Liquid and solids acceleration losses may be assumed 
negligible as a 3 ft. acceleration section was provided 
before the measuring section. Distances required for 
solids of the size and density used in this work to attain 
their constant settling velocity in water are of the order 
of 1 Hence the distance allowed for acceleration was 
more ak ample. 

To determine if a significant effect from solids-wall 
friction existed, values of the corrected pressure drop 
\P.. were plotted vs. the corresponding values of the 
static pressure drop AP,. The results are shown in Figure 
13 for the No. 107 spheres. Linear regression yields the 
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TABLE 7 
LINGAR REGRESSION OF CORRECTED PRESSURE Drop AP, vs. 
Static PRESSURE Drop AP, FoR COCURRENT COUNTERGRAVITY 
Fiow or No. 107 anp No. 113 GLAss SPHERES AND WATER 


AP. = c AP, 
No. 107 Spheres No. 113 Spheres 

D, = 0.0184 in. D, = 0.00396 in. 

pd = 155.0 1b./cu.ft. pd = 153.1 lb./cu.ft. 
c 1.030 1.008 
Se 0.02664 0.01152 
Degrees of Freedom 43 27 

© — 1,000 
t. = _ 1.11 0.703 
Sb 

trom, from tables (13) 1.68 1.70 


straight line through the origin shown on the graph. The 
slope c of the line, the standard deviation s, of the slope, 
and value of the “Student” t calculated from the equation: 


c — 1.000 
“Student” t = t, = Dwell creae eee 


Se 


0 


are also listed. Values are also shown for the No. 113 
spheres. For both systems, the hypothesis c = 1.000 can 
not be rejected at the 10% risk level (a = 10%) as the 
values of ta calculated from equation 6 are less than t,, 
Thus, it is concluded that the corrected pressure drop 
AP. equals the static pressure drop AP, for cocurrent 
countergravity flow of the two systems studied. In these 
experiments no significant contribution to the total pres- 
sure drop was therefore made by solids-wall friction. 


are given in Table 7. The corresponding values of t,, 


THE PACKED BED AS AN UPPER LIMIT FOR 
COCURRENT COUNTERGRAVITY FLOW 


The holdups (1—e), corresponding to the packed 
beds produced by slowly decrea asing the liquid flow rates 
in the batch fluidization experiments are shown in 
Table 2. A holdup exceeding the fully packed bed 
(1—e), can not exist in a uniformly fluidized system. 
Consequently, the packed bed constitutes an upper limit 
for the holdup of a fluidized system. Since the theory 
assumes uniform fluidization, the magnitude of the hold- 
up which it predicts is equal to or less than (1—e),. 
Thus a holdup equal to (1 — €)p is the upper limit lon 
ccourrent countergravity flow in a uniformly fluidized 
system. Values exceeding this can only be realized in a 
heterogeneously fluidized system or w ‘hen the solids are 
in plug flow. 

The theory says nothing about partially fluidized 
systems. Completely and partially fluidized systems may 
be distinguished according to the degree of particle- 
particle contact. In completely fluidized systems the only 
particle-particle contact is that which occurs when two 
particles collide as a result of their individual motions 
within the fluid. In partially fluidized systems particle- 
particle contact also occurs in slugs or agglomer rations of 
many particles where the duration of contact between 
two particles is considerably greater than the collision 
contact time. In partially fluidized systems, a particle may 
be mechanically supported by one or more other particles 
in contrast to support solely by virtue of the fluid flow 
around the particle. Thus, mechanical pressure applied by 
the fluid to an agglomerate in cocurrent countergravity 
flow would tend to compact the particles in the ag- 
glomerate. It is therefore possible that, due to this com- 
paction, holdups greater than (1—e), can result in 
partially fluidized systems. 

Several of experimental runs made with the No. 113 
spheres resulted in values of holdup greater than (1 — e),. 
For values of holdup less than (1—.e), random motion 
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friction vs. static pressure drop for cocurrent counter- 
gravity flow of No. 107 glass spheres and water. 


of the individual particles were observed. This random 
motion, superimposed on the net flow of particles up 
the column, is characteristic of complete fluidization. For 
three runs in which values of holdup greater than 
(1—e), were observed the random motion of the par- 
ticles ceased and packed masses of particles (agglom- 
erates) flowed intermittently as a unit up the column. 
The resulting condition resembled plug flow in which 
bubbles of liquid separated plugs of packed particles. 

These results support the requirement of the theory 
that the existence of a packed bed is the upper limit for 
the cocurrent countergravity flow in a completely fluid- 
ized system. 


Significance of the fluidization velocity 

In previous usage it has been commonly assumed that 
a unique fluidization velocity exists, characteristic only 
of the given fluid-solid system. As defined, it has been 
measured as the value of the superficial fluid velocity at 
the holdup of a fully packed bed (1— e), in a mechani- 
cally supported fluidized system (batch fluidization) with 
no net flow of particles through the fluidization zone. 

However, the present theory and experimental results 
demonstrate that a spectrum of fluidization velocities 
exists for any fluid-solid system. The fluidization velocity, 
corresponding to the initiation of fluidization depends 
upon the particle velocity V4’. As oper ation in the bot- 
tom restrained cocurrent countergravity flow region cor- 
responds to the behavior of a bottom fed fluidized bed 
which overflows the top of the column, the fluidization 
velocity may vary between two limits; at one extreme, 
that value of the superficial fluid velocity V,’ at the hold- 
up of the packed bed corresponding to ween net particle 
velocity and at the other extreme that of the single 
particle transport velocity for the critical particle velo- 
city. The fluidization velocity for any net particle flow 
between zero and the critical value is the value of the 
superficial fluid velocity at the intersection of the con- 
stant particle velocity line and the packed bed line on 
the generalized operating diagram. 

Thus the theory suggests a redefinition of the term 
“fluidization velocity” to allow for the dependence of 
the fluidization velocity on particle flow rate in bottom 
restrained cocurrent countergravity flow systems. 
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Nomenclature 


constants in equation 3 

slope of AP. vs. AP; line 

total pressure drop, lb./sq.ft. 

corrected pressure drop, lb./sq.ft. 

liquid-wall friction pressure drop, lb/sq.ft. 

static pressure drop, lb./sq.ft. 

estimate of standard deviation of bi, bs, and c 
c — 1.000 

“Student” t = or 


. 2, == 
Sc (s beocur 


(Beocur so Dpatch ) 


S*bbatch nda 
continuous phase (fluid) velocity, ft./min. 
superficial continuous phase velocity, ft./min. 
dispersed phase (particle) velocity, ft./min. 
superficial dispersed phase velocity, ft./min. 

slip velocity, ft./min. 


distance over which pressure drops are measured, ft. 


% risk 

void fraction 

holdup of fluidized bed 

holdup of packed bed 

density of continuous phase, Ib./cu.ft. 
density of dispersed phase, Ib/cu.ft. 
viscosity of continuous phase, centipoises 
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Air Agitation and Pachuca Tanks 


A. G. W. 


An expression for the energy transferred from 
air bubbling through the contents of a tank is de- 
veloped as a function of the fluid depth in the tank, 
the specific gravity of the tank contents, the vapor 
pressure of the solution, atmospheric pressure, and 
the volume of air used. The variation of the energy 
transfer with height from the bottom of the tank is 
also developed. 


It is shown that a portion of the energy trans- 
ferred is used up in small-scale turbulence and eddy 
effects, and the remainder is available for genera- 
tion of a circulation pattern in the contents of the 
tank. 


The information developed is used to determine 
the pulp circulation pattern in pachuca tanks of 
various design. The functions fulfilled by pachuca 
tanks are discussed, and the considerations affecting 
a choice of tank described. 


in currents for the agitation of the fluid contents 
As a tank have long been used in industry, some- 
times as a temporary measure, sometimes as a permanent 
one. The satisfaction achieved with air as an agitating 
medium varies from one application to another. C. S. 
Quillen “) says that “basically, air agitation is inefficient, 
yet it fills a real need in applications such as corrosive 
liquid jobs, fluidizing high-percent-solids slurries, or agi- 
tating shallow vessels where mechanical mixing is im- 
practical”. Another interesting new use for air agitation 
has developed in the past few years, in the aeration of 
harbors and channels in the winter time for the pre- 
vention of ice. 


Some decades ago a special application of air agitation 
was invented in the metallurgical industry, namely, the 
Brown, or pachuca tank. The pachuca tank is a vessel 
with a conical bottom (usually with a 60° included 
angle) having air introduced at the apex of the cone. 
The conical bottom is intended to direct settled solids 
into a region from which they may be returned to the 
top of the tank by the fluid currents developed by the 
air. This type of tank was first used at Komata Reefs in 
New Zealand. A partial record of its development is 
given by von Bernewitz ). The leaching of gold ores in 
pachuca tanks became common in South Africa early in 
this century. When the uranium industry was established 


iManuscript received February 24, 1958. 
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in that country in this decade, the leaching vessel every- 
where chosen was the pachuca tank. 


In Canada, the uranium industry has recently found 
an interest in the pachuca tank for its ore leaching op- 
erations. The first Canadian mill to use this type of tank 
was the Beaverlodge mill of Eldorado Mining and Re- 
fining Limited. At Beaverlodge, uranium ore is leached 
in alkaline solution, and the oxidant requirements of the 
process are supplied by aeration. In Canadian uranium 
operations employing ‘the acid leach process, Bicroft 
Uranium Mines Limited was the first to make use of the 
pachuca tank. Some of the Blind River mills are also 
using it. In the acid leach process, there is no process 
benefit to be derived from aeration; the whole purpose 
of the air is for agitation and suspension of solids. For 
this purpose, the pachuca tank has the attraction of ex- 
posing no moving or mechanical parts to the acid pulp. 


Despite the widespread use of air agitation, there is 
very little literature of an analytical nature available on 
the subject. Kaufman), in a few lines, has suggested 
practical amounts of air for various intensities of agita- 
tion, referring specifically to the petroleum industry. 
King *) is concerned mostly with reporting actual prac- 
tice in pachuca tanks. 


With the revival of Canadian interest in the pachuca 
tank, some analysis of the principles of air agitation and 
of the nature of pachuca tank agitation in particular 
seems desirable. This paper presents such an analysis on 
the basis of the energy transfer in the tank. The argu- 
ments, discussions, and conclusions presented are almost 
entirely rationalizations from principles and, where neces- 
sary, assumptions. The development of experimental 
verification of the numerical arguments is highly desir- 
able. Free-airlift circulation (as defined in following sec- 
tions) is an unexplored field for experimental study, and 
it is hoped that the discussion here will stimulate interest 
in it. 


Energy transfer in air agitation 


A given volume of air under pressure contains an 
amount of energy from which useful work may be ob- 
tained under certain conditions. The work is derived 
from the conversion of some of the heat energy of the 
gas (adiabatic expansion) or of the surroundings (iso- 
thermal expansion) to mechanical form. The process of 
conversion is called here “energy transfer” from the air, 
and the term is taken to mean the conversion of the 
potential energy of the compressed gas to other forms 
of energy. 
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The agitation developed in an air-agitated vessel re- 
sults from the transfer of energy from the air to the tank 
contents as the air expands in rising through the vessel. 
In analyzing the source of the agitation energy, it is 
helpful to start with the simplest case — that in which 
the solution vapor pressure is low enough that it may be 
neglected. The method will then be broadened to include 
vapor pressure contributions. Two assumptions are made 
in the following discussion. The first is that the expan- 
sion of the air in a tankful of fluid is isothermal. This 
assumption is probably very close to being a statement 
of fact, and actual minor deviations from it will affect 
the argument and conclusions little. The second is that 
the analysis may reasonably start w ith air at the tank 
bottom. Nozzle and jet effects, and the kinetic energy 
of the entering air, are ignored. 


In the simplest case, the analysis may be based « 
study of the air history in the tank or on that of Pn 
pulp. Both approaches are used below. 


Method A. A given volume of free air exhausted at 
the top of the tank has transferred energy in expanding 
through the pulp equal to the isothermal work of com- 
pression to the pressure at the bottom of the tank. The 
work equation is 


W = 144 PIV: loge 5 ++ Salita 
1 
where W is the rate of energy transfer from air to pulp, 
in ft.lb/min 

P is pressure, in lb./sq.in. (absolute) 

V is air volume rate in cu.ft. per min. (unless 
stated differently air volumes in this paper 
refer to air at barometric pressure and at the 
temperature of the tank contents) 

subscript 1 refers to the tank top. 
subscript 2 refers to the tank bottom. 


The work equation is more conveniently written 


W = 144 P,V, loge ee ee poh eee eee 
‘ 2.31 P, 
where s is the specific gravity 
D is the depth at which air is released into the 

tank in ft. 
As an example of the application of the work equation, 
consider a tank containing 45 ft. of pulp at sp.gr. 1.6. 
For each 100 cu.ft. per minute of free air at 14.7 Ibs/sq. 
in. (absolute), the rate of energy transfer from the air 
amounts to 242,000 ft.lb./min. or 7.34 H.P. 


Method B. The rise of a volume of air from the 
bottom of the tank to the top is, conceptually, equivalent 
to the fall of an equal volume of pulp from the top of 
the tank to the bottom. Let L Ibs. be the weight of pulp 
of volume V. Then, L = 62.4 s V, and the energy con- 
verted in the fall from potential to mechanical form is 
w= H where H is the height from the tank bot- 
tom, in ft. Then, dW = 62.4 s V,, dH where subscript H 
refers to conditions at height H ft. Assuming the perfect 
gas law, and substituting for V in the above equation, 


dW 144 PV, s 


= (3) 
dH 2.31 Pi +s (D — H) 


This equation gives the amount of energy transfer per 
foot of depth at any depth. It is a simplified form of the 
basic air agitation work equation developed below. In- 
tegration of the differential equation(3) gives the same 
overall work equation as given in Method A. The two 
methods are thus equivalent in their answers. 
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The general air agitation work equation 


In the case where solution vapor pressure is appre- 
ciable, the volume of gas passing through the pulp is 
increased because of evaporation into the gas stream, 
and the effect of the increased volume must be token 
into account in the work equation. The events taking 
place in the tank may be described qualitatively as fol- 
lows: 


1. Air at the tank bottom is initially at a pressure 
corresponding with the static head of pulp above it. 
2. Liquid evaporates into the air, decreasing the par- 
tial pressure of the air until saturation is reached. The 
heat of evaporation is assumed to be replaced from 
an external source, and the departure from isothermal 
conditions at any point in the tank is taken to be 
minor. 

3. The air and vapor mixture rises through the tank, 
increasing amounts of vapor being carried along to 
maintain saturation. 

4. At exhaust, the partial pressure of the air is less 
than atmospheric pressure. 


It is assumed that equilibrium saturation exists at all 
depths, and that perfect gas laws apply. 

The work equation for the rate of energy transfer 
from air to pulp per foot of tank depth at any depth is 


dW f PiVo 
= 1445 
dH {2.31 (Pi — p) + s(D — H){ 
f 2:31 p 
yi + (4) 


| 2.31(P:— p) +s 


where p is the solution vapor pressure, lb./sq. in. 
subscript o refers to dry air at tank top pressure 
and at tank temperature. 


aa 


This is the general air agitation work equation. Integra- 
tion will give the overall total rate of energy transfer 
from air in the tank. If the vapor pressure p is zero, the 
equation simplifies to the one developed earlier for that 
case, V, then being equal to V, 


Energetic characteristics of air agitation 

In a mechanically-agitated tank, the mechanical energy 
of the mechanism is transferred from the impeller to the 
tank contents in a zone restricted in the vertical plane 
to the vicinity of the impeller. The air- agitated vessel 
differs profoundly in having the energy input diffused 
throughout the depth of the tank. Figure 1 presents some 
differential energy transfer rates calculated from Eq. 3 
and Eq. 4 for various conditions. 

Three air-agitation conditions are shown in the figure 
along with a typical curve for mechanical agitation. With 
the mechanical impeller the greatest agitation intensity 
is near the bottom, presumably where greatest intensity 
is needed for sand suspension. For air agitation, however, 
the intensity is lowest at the bottom, and can increase 
greatly with approach to the fluid surface. Particularly 
if the fluid vapor pressure is appreciable, the dispropor- 
tionate energy transfer rate toward the tank top becomes 
very marked. A given volume of dry free air is more 
active throughout the tank when assisted by a high vapor 
pressure. Conditions at the top of the vessel, however, 
can become explosively violent. 

When the top of an air agitated vessel is under pres- 
sure, the energy potential carried into the bottom of the 
tank in the air is partly transferred to the tank contents 
and partly carried out in the pressurized exhaust. The 
amount of free air required for a given energy transfer 
rate to the pulp is thus increased. Unless air or gas is 
necessary to the process, therefore, air agitation must 
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become economically less and less attractive as the top 
pressure increases. 

What becomes of the energy transferred from air 
to pulp depends greatly on tank and air inlet geometry. 
In a mechanically agitated tank (suppose a vertical cylin- 
der with a turbine impeller), the energy history is usually 
completely changed by the addition of baffles on the tank 
wall. Without baffles, most of the energy is dissipated 
in eddy turbulence at the tank walls because of the swirl 
of the fluid contents. With baffles, most of the energy 
appears first in fluid currents from the impeller, which 
then are degraded to eddy currents and heat throughout 
the tank. An analogous situation exists with air agitation, 
gross fluid currents or generalized eddy turbulence may, 
develop depending on the gas distribution over the tank 
bottom. A local gas addition will result in gross current 
analogous to those in the baffled, mechanically-agitated 
tank. 


Loss of energy to turbulence and eddying 


As gas bubbles pass through a tank, a portion of their 
energy potential is lost because of bubble slip. As the 
bubbles rise, fluid is displaced from ahead of them to 
behind them, with resulting small-scale turbulence and 
eddying. Depending upon the geometric relationships of 
the tank and the air current, another portion of the 
energy potential may be converted to gross fluid flow as 
discussed above. The proportionate amounts of the 
energy potential that are lost in turbulence and that pro- 
duce gross flow can readily be determined. 


Let the relative velocity of bubbles and pulp be 
z ft./sec. and the velocity of the pulp stream in which 
the bubbles move be y ft./sec. The time required for = 
bubbles to rise through a height dH ft. is dH/(z+ y 
seconds, and during this time the slip or the excess x. 
tance covered by the bubbles as compared with the pulp 
. dH oa 
is ay 2 ft. Using the argument of Method B above, the 
zt+y = 
loss of energy to turbulence and eddying is then 


dW1 dW 1 


= ae ’ (5) 


dH dH (1+ y z) 


where W, is the turbulent energy loss caused by bubble 
slip, in ft.lb./min. : , 

The above equation shows that a certain fraction of 
the potential energy of the air is lost to small-scale tur- 
bulence, and the fraction is 1/(1+ y/z) of the total 
energy transfer. If there is no gross pulp flow (y=0), 
the fraction becomes unity, and the whole of the energy 
potential is converted to minor eddy currents. This is, 
for example, what occurs where the air is distributed 
uniformly across the whole tank cross section. The 
higher the pulp flow rate, the lower is the fraction of the 
energy lost to slip. The other factor affecting the fraction 
is the rate of bubble rise. 

Taggart ©) reports that 2 mm. air bubbles rise through 
a pulp at a rate of 150 mm./sec. or about 6 in./sec. U.S. 
Patent 2,740,707 ‘“) quotes the rate of rise of gas bubbles 
in a specialized type of leaching tower as 4 to 8 in./sec. 
Using a value for z of 0.5 ft./sec. in Eq. 5, the fraction 
of the energy transfer lost to bubble slip is 1/(1 + 2y). 

It is worth noting that the loss fraction is about 15% 
in pachuca tanks with no centre column, and that it could 
probably be greatly reduced by providing finer bubbles. 
Leibson et al), working with orifices up to 4” in dia- 
meter, have shown that mean bubble diameter is a func- 
tion of the Reynolds number in the orifice. E xperimental 
study of the methods for introducing air into pachuca 
tanks might result in appreciable gains in efficiency. 
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Figure 1 — Variation of energy transfer rate throughout 


a tank. 

Line A—Open top, 14.7 Ib./sq.in. barometer, vapour pres- 
sure zero, 100 cu.ft./min. free air. 

Line B—Top pressure 215 Ib./sq.in. (absolute), vapour 
pressure zero, 850 cu.ft./min. free air. 

Line C—Open top, 14.7 lb./sq.in. barometer, vapour pres- 
sure 10 Ib./sq.in., 30 cu.ft./min. free air. 

Line D—Mechanical turbine, 1 ft. deep blades. 

Notes — Integrated areas below all curves are equal. 
Specific gravity 1.6 Pulp depth 45 ft. 


Pulp circulation in pachuca tanks 

Using the principles and equations developed above, 
it is possible to estimate the nature of the pulp circulation 
in pachuca tanks, and to compare the effectiveness and 
economies of various types. In all the e: xamples, a vessel 
having a 45 ft. depth of pulp of specific gravity 1.6 is 
considered. A standard pachuca tank would be 22’6” i 
diameter, with a 60° cone at the bottom, and would nor- 
mally use about 300 cu.ft./min. of free air“). Three types 
of tank have evolved since the original installation at 
Komata Reefs; each of these, and a ‘proposed new one, 
is discussed. The four types are sketched in Figure 2. 

The original pachuca tank had in it a central column 
for almost the full tank depth (A, Figure 2). Air is intro- 
duced at the bottom of the column, causing pulp to rise 
in the column, so developing circulation in the tank. The 
second type of tank was evolved by a arg in 
cyanide practice, as recorded by von Bernewitz ‘*). When 
the upper sections of the full- length airlift coliam were 
successively raised slightly and then removed, improved 
extraction of gold was obtained. This discov ery resulted 
in the development of a tank having the centre column 
extending as a stub only to the top of the conical section 
of the tank (B, Figure 2). The third type of tank in 
actual operation is that used almost universally in the 
South African uranium plants. In it, the centre column 
has been a removed (C, Figure 2). A fourth 
type of vessel proposed here has a full centre column 
but differs from the Komata Reefs tank in the point of 
air admission (D, Figure 2). A consideration of the func- 
tions required of the processes going on in a pachuca 
tank is useful in the analysis of performance. 
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Figure 2 — Sketch showing types of a pachuca tank 
discussed. 


cscs # tame 


A—Full-centre-column, or Komata Reefs tank. 


The first requirement of a pachuca tank is that it 
maintain the suspension of solids in liquid that is fed it 
— ie. that it not allow sand to settle cumulatively. The 
second is that it may be required to aerate the pulp con- 
tained in it. Finally, it may be required to scrub. films 
off the solid particles present. The four types of tank 
described may fulfill each of these functions differently, 
and one object of the analysis below is to examine these 
matters, 


A. The Komata Reefs, or full-centre-column tank. 
In the machine with a full centre column, the energy 
transferred from the gas is used to produce entrance, 
friction, and velocity heads associated with pulp flow in 
the column. Each of these hydraulic heads is a function 
of the fluid velocity, y, with the following values “): 


Entrance head, for a Borda-type 
/> 


entry 0.83 y*/2g 
Friction head, for a friction factor 

of 0.03, a tube length of 45 ft., 45 

and a tube diameter of 1.5 ft. 0.03 3 y*/2¢ 

iG ae 

Velocity head 1.0 y*/2g 
Total hydraulic head, the sum of 

entrance, velocity and friction 

heads 2.76 y*/2g 


The pumping rate through a centre column may 
be easily calculated as 1180 cu, ft. per min. The calcu- 
lation is based on the energy Ben rate (242,000 ft. 
lb./min. per 100 cu. ft./min. of air passing through 45 ft. 
of pulp of specific gravity 1.6), the relationship between 
pulp velocity and pumping rate, and bubble slip allow- 
ance from Eq. 5 5. 

Cook and Waters“) have reported extensive tests on 
pumping rates through submerged gas-lift circulators 
used in tanks for the storage of radioactive wastes. They 
have also proposed a procedure for the calculation of the 
flow through such circulators, using in the procedure 
certain constants derived from their experimental work. 
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B—Stub-column tank. 
C—Free-airlift tank. 
D—Full centre column, with shallow air introduction. 


Their method predicts a flow of about 960 cu. ft. per 
min. through the 18 in. tube in the tank being considered 
here. Inasmuch as Cook and Waters feel their calculation 
method to be good to only + 25%, the agreement with 
1180 cu. ft. per min. above is not unreasonable. The 
principle reason for the difference is that where 2.73 
velocity head is used for the pumping head above, Cook 
and Waters suggest about 4.6 velocity head. Their 4.6 
velocity head is based on tests that involved an air pipe 
passing down through the length of the column, and on 
a gas distributor more restrictive of pulp flow than nor- 
mal in pachuca tank practice. Cook and Waters method 
predicts a flow of 1180 cu. ft. per min. with 300 cu. ft. 
per min. of air if the tube diameter is taken a few inches 
over 18 1n. 

The flow of 1180 cu. ft. per min. is plotted in 
Figure 3 as line A. 


B. The stub-column tank. In the tank having a 
centre column extending only to the top of the cone, 
pulp energetics similar to those in a full column are 
experienced in the stub. Above it a free airlift exists. 
For a 20 ft. column of 18 in. diameter, the pumping rate 
in the column is calculated to be 800 cu. ft. per min., 
compared with 1180 cu. ft. per min. through a full 
column. 


The behaviour of the stub-column tank above the 
column is similar to that of the free-airlift tank. Since 
the velocity of discharge of pulp from the column is 
close to that which is suggested in the following section 
for the central axis of a free airlift, it is assumed that 
the pumping curve above the column parallels that of 
the free-airlift tank. The flow is plotted in this way as 
line B in Figure 3. 


The free-airlift tank. Free-airlift tanks in. op- 
eration with an air rate of 300 cu. ft. per min. are ob- 
served to have the air released in an area of 20 to 60 
sq. ft. The release area is raised about 12 in. above the 
level of other parts of the pulp surface. Pulp which has 
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Figure 3 — Pumping rates in various types of pachuca 
tank. 


Line A—Komata Reefs tank, 18 in. X 45 ft. column. 
Line B—Stub-column tank, 18 in. X 20 ft. column. 
Line C—Free-airlift tank, no column. 


Notes—Air rate 300 cu.ft./min. Specific gravity 1.6. Pulp 
depth 45 ft. 
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been raised is observed to flow across the top and to the 
sides of the release area. Finally, release is not a steady 
phenomenon, but, rather, the release area “wanders” at 
random around the top pulp surface. These facts repre- 
sent the extent of the experimental data available for 
the estimation of pulp circulation rate in the free-airlift 
tank. Using one of the facts (the 12 in. elevation of the 
release area), the other two (the extent of the release 
area, and its wandering) are confirmed in the analysis 
below. Apart from these observed facts, the analysis is 
entirely a rational one; the various assumptions that have 
been necessary are discussed explicitly. As noted in the 
introduction, experimental study of these assumptions 
and of the mechanisms proposed is highly desirable. 
Meanwhile, the circulation pattern developed in the ra- 
tional analysis is shown in Figure 4. 

Where the upflowing pulp springs higher than the 
general pulp level in the tank, and is not hindered or 
restricted by cross- flowing pulp, the pumping head is 
equal to the height of spring. This statement is equival- 
ent to saying that, in a cross section at the level of the 
static general pulp — air interface, the pressure below the 
free-springing area (or point) is atmospheric. At other 
points, the upflowing pulp is hindered by cross-flowing 
pulp, and the pumping head at these points is the sum 
of the rising velocity head and the depth of cross-flowing 
pulp. 

By rough analogy with work on free jets (for ex- 
ample, Folsom and Ferguson (10) ), it is assumed that the 
velocity distribution around an instantaneous central axis 
of the rising stream is such that « = 1 — r where uw is the 
velocity on the given axis over the velocity on the central 
axis, 1 is the full radius of the rising column, and r is the 
distance of the given axis from the central axis. Such a 
distribution of velocities can be calculated to imply an 
average velocity equal to one-third that on the central 
axis. 

These concepts are sketched in Figure 5, using flow 
rates and dimensions from the following assumptions and 
calculations. The sketch shows a condition that must 
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Figure 4 — Representation of flow in a free-airlift 
pachuca tank. 


Tank size 22 ft. 6 in. X 45 ft. of pulp. Air rate 300 cu.ft./ 
min. Specific gravity 1.6. Solution vapour pressure zero. 


Number of flow lines shown entering the rising pulp 
stream per foot of height is proportional to the volume 
rate of flow developed from each level. 


be unstable because of the “pressure blister” above the 
rising pulp stream. This pressure blister results from the 
presence of cross flowing pulp at the surface, and is the 
cause of the observed wandering of the release area over 
the tank surface. The sketch shows also that pumping 
head is made up of a combination of velocity head and 
pressure head. 

As a working approximation, it is proposed to assume 
that the pumping head is 12 in. across the whole rising 
pulp stream. (At the central axis, this head has only a 
velocity component; at the periphery of the rising stream, 
the 12 in. head is entirely pressure head resulting from 
the cross flowing pulp stream.) 
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Figure 5 — Sketch of relationships in the release area of 
a free airlift. 


With a pumping head established, the question arises 
of what portion of the energy transferred from the air 
results in “useful” pulp flow. Useful flow is considered 
to be that in a vertical direction. Loss of energy to non- 
useful pulp movement results from two sources; bubble 
slip produces small-scale turbulence as already discussed, 
and horizontal movement of pulp from the surroundings 
into the rising stream produces large-scale turbulence. 
The latter can be approximated from a consideration of 
the velocity head involved in horizontal pulp flow. Using 
figures which are of the correct order of magnitude, to 
pass 6,000 cu.ft./min. of pulp horizontally and uniformly 
into an inverted cone of height 45 ft. and base radius 
7 ft. requires an entry velocity of less than 0.2 ft./sec. 
This represents a velocity head of less than 1% of the 
total vertical pumping head and so may be properly 
ignored. 

The last condition required to define the flow regime 
is the vertical velocity distribution. If it is assumed that 
the mechanism of the transfer of energy from air to 
pulp is similar at all levels in the tank and that terminal 
bubble velocities are established immediately at the tank 
bottom, then the assumption of uniform upward pulp 
velocity on any given axis is reasonable. It follows that 
pulp flow at any level is proportionate to the energy 
transfer to that level. 

For a velocity head of 1.0 ft. on the central axis, the 
corresponding velocity is 8 ft./sec. and the average up- 
flow velocity one- -third of this or 2.67 ft./sec. For a rela- 
tive bubble velocity z of 0.5 ft./sec. it has been verified 
that a refined calculation of slip loss based on velocity 
distributions gives a figure very close to that based on the 
average velocity. Using the latter method, slip loss is 


0.5 , ; : 
100% or 16% of the available power. For 








an energy transfer rate of 726,000 ft.lb/min. per 300 
cu.ft./min. of air, the full flow developed in the airlift 
726000 0.84 


1.0 
per min. The area of the rising stream at the top of the 
tank is then 38 sq. ft., and the diameter 7 ft. The flow at 
various levels is plotted in Figure 3 as line C. Figure 4 
is a scale representation of a cross section through a free 
airlift tank, based on the details of the analysis above. 


is then = 610,000 Ib. per min. or 6100 cu.ft. 


D. A proposed new pachuea tank. In the conven- 
tional pachuca tank, to the author’s knowledge, air is 
always introduced at the bottom. (The tanks at Eldor- 
ado’s Beaverlodge operation have a split air addition, 
part going in at the bottom of a centre column and part 
into other levels of the tank to produce aeration.) A 
consideration of the methods used in this paper suggests 
that air addition need not be made at the bottom. In 
some cases, the ability to add air at other than the tank 
bottom may have ‘distinct advantages. For ex ‘ample, 
where tanks of various heights are used, power economies 
may be had if air need not be compressed to a pressure 
suitable for the deepest tank. As an example of the 
potential of the new tank proposed, 300 cu.ft./min. of 
air introduced into a full centre column 35 ft. above the 
bottom, in a tank otherwise like the Komata Reefs tank 
discussed above, would develop a flow of 910 cu.ft./min. 
The air pressure at point of injection would be only 

lb./sq.in. (gauge) compared with 31 Ib./sq.in. (gauge) 
for conventional practice. 


Choice of type of tank 
As noted in the previous section, a vessel used for 
leaching ores may be required to perform one or more 
of several functions. Each of these functions is affected 
by different factors as follows: 
a. Maintenance of suspension of solids requires either 
a positive flow from some low level in the tank or a 
general random turbulence in the tank. If a positive 
current is depended upon for suspension, it must 
carry up all the sand that falls into its source. In most 
cases, though not universally, this one mechanical 
function at least is required of a pachuca tank. 
b. Aeration of the pulp may be a composite function 
of the volume of pulp in intimate contact with air, 
the volume of air in intimate contact with pulp, and 
the pressure. Aeration is important to the cyanidation 
of gold ores, and played an important role in the de- 
velopment of the pachuca tank in South Africa. 


c. Scrubbing off of films from the solid particles is 
probably most affected by the rate of shear associated 
with agitation. Some ores may require no scrubbing 
action, and others may require some minimum amount 
to remove, for example, gypsum films, It has some- 
times been suggested that liquid films are important 
to ore leaching rates, though the author disagrees with 
this contention, at least as applied to the normal leach- 
ing of uranium ores, 


In the metallurgical industry at present, there is far 
less than enough data to allow a quantitative matching 
of the characteristics of the various pachuca tanks with 
the needs of a particular ore leaching operation. A quali- 
tative assessment of the characteristics of each type of 
tank will, however, assist in general considerations of 
their applic: ability and perhaps in attaching meaning to 
comparative leaching tests in various kinds of vessel. 

For applications requiring particle scrubbing, the 
mechanical agitator should be expected to be far superior 
to any air agitated tank. Scrubbing is related to the 
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amount of shear, which in turn is a function of relative 
velocities and of energy transfer rates. Of the pachuca 
tanks, the Komata Reefs tank develops the highest rela- 
tive velocities, though they are far less than can exist in 
mechanically-agitated tanks. This tank may be better 
than the others for light particle scrubbing. 

Aeration is obviously best in the free-airlift and sub- 
column tanks. The improved gold extraction recorded by 
von Bernewitz‘*) in the progressive conversion of a 
Komata Reefs to a stub-column tank was _ probably 
caused by improved aeration. 

For sand suspension, the Komata Reefs and _ stub- 
column tanks would appear to be superior to the free 
airlift tank. Figure 3 shows that the pumping rate in the 
free airlift does not equal that in a stub column until a 
height of about 10 ft. above the bottom is reached. The 
free airlift tank, therefore, will have a tendency to ac- 
cumulate sand at the bottom to a greater extent ‘than the 
tanks having centre columns. However, for pulps that 
sand out severely and rapidly, the free airlift tank has 
a natural equilibrium which may be important. So long 
as the pumping rate through a full centre column is 
enough to raise all the sand that settles to the bottom, 
the bottom will stay clean. If the pumping rate is not 
enough to raise it all, a continuing sand build-up will 
occur until the tank is full of sand or sand is drained 
from the lower levels. The free airlift tank on the other 
hand will have a sandy bottom of depth varying with 
sanding conditions. The stub-column tank has partly the 
nature of both of the others, unstable to the top of the 
column and with varying equilibrium above it. For most 
pulps handled in pachuca tanks rapid and severe sand 
settlement (rapid in relation to pumping rates) is not a 
problem. The Komata Reefs and stub-column tanks then 
have the advantage of keeping the bottom quite clean 
with less air than the free airlift tank will do. 


Scale-up and analogy 

The usual pattern of testing for ore leaching opera- 
tions ignores possible scale-up problems so far as leaching 
is concerned. Bench tests may be done in a beaker with 
a laboratory propeller mixer of, usually, undefined char- 
acteristics. Pilot plant work may involve a 45 gallon drum 
and } h.p. portable mixer The full scale plant may use 
pachuca tanks of any design, or mechanical agitators, or 
air lifts combined with sw eep agitators (Dorr). One stage 
is seldom related to the other insofar as mixing intensity 
is concerned. It is fortunate that, at least in most uranium 
and gold ope rations, the process is not particularly sensi- 
tive to mixing intensity. Some of the laboratory-plant 
differences which are actually experienced, however, may 
be caused partly by such a factor. It should at least be 
thought of in the approach to a new ore or new process. 


In comparing pachuca tank operations from place to 
place, and in designing new installations, two different 
dimensional factors are commonly used. One refers to 
air consumption in terms of cu.ft./min. per ton of solids 
contained in the tank. This is equivalent, for standardized 
pulp and solids density conditions, to using cu.ft./min. of 
air per cu. ft. of tank volume. The other basis refers to 
air in terms of cu.ft./min. per sq. ft. of tank top area. The 
first implies the air requirements vary as the area and 
the depth; the second implies that tank depth is unim- 
portant, 

Figure 6 shows differential energy transfer rates as a 
function of the depth of tank. If it is assumed (as in the 
previous section) that one unit of power will produce 
one unit of flow, then Figure 6 can be used directly for 
guidance in scale-up. For equal agitation at the tank tops, 
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Figure 6 — Variation of energy transfer rates with 


tank depth. 
Line A—Differential transfer rate at the bottom of the 
tank. 
Line B—Overall average transfer rate. 
Line C—Differential transfer rate at the top of the tank. 


Notes—Air rate is 100 cu.ft./min. in each case. Specific 
gravity 1.6. Solution vapour pressure zero. 


tank depth is not a factor; referring to air consumption 
in terms of cu.ft./min. per sq.ft. of tank top area may 
then, properly, be used. For equal agitation at the tank 
bottoms, a shallower tank requires less air than a deep 
tank, but the relationship is not an inverse proportional 
one. It is, in fact, such that 
Va Ve 
= z (6) 
34 + s. Da 34 + s. Dg 


where subscripts A and B refer to tanks A & B. 
Whether to refer air consumption to the volume of the 
tank or to the tank top area depends, then, on the con- 
ditions to be defined in the tank. Since air-agitated tanks 
in the met allurgical industry are usually used primarily 
for solid suspension, bottom conditions are most im- 
portant. Perfect scale-up and analogy in that case exists 
where the cu.ft./min. per sq. ft. per ft. absolute static 
pressure at tank bottom is maintained equal from tank 
to tank. Small-scale testing of pachuca tanks can be made 
with confidence on such a basis. When solution vapor 
pressures are appreciable, they must be allowed for in 
accord with the general air agitation work equation. 

Kaufman (*) has given the following quantities of air 
per sq.ft. of tank cross-sectional area for a liquid depth 
of 9 ft.: 

Cu.ft./min. free air 


Moderate agitation 0.65 
Complete agitation 1.3 
Violent agitation 3.1 


He also says that a 3 ft. depth of solution requires tw ice 
the quantity of air “for the same degree of agitation”. 
Since Kaufman presents no foundation for his statement, 
it is assumed that his figures are based on the appearance 
of the surface of the fluid. Surface appearance is a direct 
function of the total fluid flow rising to the surface; for 
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equal pumping heads, it is then a function of the total 
energy transfer rate in the tank, Eq. 2 must then have 
equal value in the two tanks. It can be calculated from 
Eq. 2 that a 3 ft. tank should require two and a half times 
as much air as a 9 ft. tank. The basis assumed above for 
Kaufman’s statement about the effect of depth is, there- 
fore, confirmed, if 24 times by calculation is accepted as 
agreeing with Kaufman’s “twice”. Kaufman’s undefined 
“degree of agitation”, then, is simply a reference to the 
appearance of the surface of the fluid. His 9 ft. and 3 ft. 
tanks should be compared on the more relevant criteria 
of amount of flow generated from any given level. The 
3 ft. tank on this basis requires the same amount of air 
as the 9 ft. tank for equal flow generated per ft. of depth 
at the top, or § as much for equal flow generated per ft. 


of depth at the bottom. 
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Mixing Studies on a 
Perforated Distillation Plate 


A. I. JOHNSON? and J. MARANGOZIS: 


In establishing a relationship between the plate 
efficiency and point efficiency for perforated plates, 
splashing of the liquid was considered as an entirely 
new mechanism for the mixing of the liquid phase 
in the direction of the flow. Actual splashing 
measurements were reported and correlated here, 
under different conditions of flow rates, weir height, 
kinematic viscosity of the liquid and some aspects 
of tray geometry. 


Plate efficiencies for the desorption of ammonia 
by air were reported for a perforated plate 36 inches 
long and 3 inches wide. Concentration gradients 
across the plate were also measured for the am- 
monia desorption case and for backward distribu- 
tion of sodium chloride by splashing. These con- 
centration gradients have been related to the splash- 
ing measurements. 


_ and plate efficiencies were defined by Murphree) 
as: 


and 


respectively, 

Lewis ) assumed no mixing whatsoever in the direc 
tion of the flow and complete mixing in the vertical 
direction and derived the following relationship between 
the above efficiencies: 


kemv (3) 


where 4 is mG/L. 

Several investigators tried to modify Eq. (3) to in 
clude the effect of mixing; Gautreaux and O'Connell ©) 
utilized the “liquid pool” concept of Kirschbaum, Oliver 
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and Watson‘ proposed a recycle model and defined a 
mixing factor C which was introduced into a modified 
Eq. (3); Anderson ©) derived an equation using an “eddy 
diffusivity” term to account for back diffusion. 
However, none of the above methods have presented 
a completely satisfactory solution to the problem in hand. 


In this investigation splashing of the liquid was con- 
sidered as the mechanism of mixing; actual measurements 
of splashing were carried out and correlated as a function 
of the hydrodynamic conditions and liquid physical 
properties. 

A visual observation of the performance of the 

bubble-plate, revealed that a thin liquid laver existed on 
the bottom of the plate, while the remainder of the 
liquid holdup was in the form of eddies and spray or 
foam; these eddies extended from the liquid layer to the 
top of the foam, they swirled and they fell back again 
on the plate; in this way they caused the liquid to 
“splash” from place to place and mix with the adjacent 
liquid upstream, or downstream. Figure 3 illustrates the 
sampling set up for measuring the amount of splashing 
discussed here. 

This splashing was considered here as the major mix- 
ing mechanism, while the mixing by eddy diffusion in 
the thin liquid layer was shown to be negligible by the 
sodium chloride experiments. 


Theoretical section 

Keeping this mechanism of mixing in mind the prob- 
lem was treated by the following mathematical analysis: 

Consider the perforated plate in Figure |. A fraction 
Q,, of the liquid flow rate L, splashes from point B to a 
distance Z, upstream and enters cross section AAL, cor- 
respondingly a fraction Q, 
downstream, while the fraction (1— Q,—Q,) 
regularly on the floor of the plate. 


splashes to some distance Zy, 
travels 
Assume now, that the , and Xp can 
be related to x, by a linear relationship, Le. 


composition x 


Xt x Zi (4a XB x + ; LB (4b) 


By a material balance around a difterential slab of unttorm 
composition vertically, there results, 


+ (1 Ox Oy) | dl (Ya ‘ dG (5S) 











A106, Yn+1 


Figure 1 — Mixing on a perforated plate by the 
mechanism of splashing. 


Now, assuming that 


ogg SORE APD es ea vice wens ones (6) 


and ya+1 = constant....... 


and substituting for 


iz 
WE ie Brn bioccciscwtciwnrsead (8) 
Zo 
nG 
Mite TE Sccdstaetacoweseneuasd (9) 
L 
and 
i. Yn — Yn+1 
Eoc = appa a (10) 
Eq. (5) could be rearranged as: 
d? (yn — Yn+1) d (Ya — Yn+1) 
(Qrzr — Qpzp) a = -_ 
dz? dz 
AEoc 
i cer MN oicosrs' sn cosade (11) 
Zo 
Define the mixing factor due to liquid splashing as: 
8 = inte = Deiei 66653 sees panes (12) 
Then 
d? (yn —_ Yn+1) d (yn — Yn4+1) 
B - o Zi 
dz? dz 
AEoc 
(Sa = Watt OOs ew ickieivende’ (13) 
Zo 


This is the general differential equation describing the 
situation when liquid mixing is present. 


It is easily seen that § is a factor which should be 
measured experimentally and represents the effect of mix- 
ing on the concentration gradient. For no splashing pB=0 
and the differential Eq. (13) gives rise to Lewis’ Eq. (3) 
where plug flow is assumed. 


The complete primitive of the second order differen- 
tial equation above is of the form: 


yn — Yori = kyeAiz + kyebrz, (14) 


where 6, and 6, are the roots of the auxiliary quadratic 
equation, and k, and k, are constants defined by the 
boundary conditions, i.e. 

for 


z= Zo, (Ya — Yati) = k, + ke. 


and for 


Z = Zo, (Yn — Yntilo = k e120 of ke92Z0 pea ieacareeaee (16) 


However, it was generally found in these experiments 
that k, was such a small number, that the contribution of 
the group {jez was always negligible, in fact of the 
order of magnitude of 10-°° to 10-*°° as compared to an 
order of 10-* of the group k,e4:z 

Therefore one could assume k, = o when the com- 
plete primitive becomes: 


ig Rg WO oa sy cu sncees (17) 
In Eq. (17) 6, is given by 
1 a 
0, = -4/ TO cane ates (18) 
268 4B? Bz. 
and 
(YnJo — Yn+1 
k, = — Ratner gaan a ase Be eee (19) 
Azo 
By definition 
n avg: ~ "n+ ) 
Rey hs oe ae re acl (20) 


(y*n)o — Yn+1) 


and (Y,)ayg could be obtained by integrating y, over the 
length of the tray, i.e. 


Zo wo 
1 1 5 
(yn)avg. = (yn) dz = — | (kye%2 + yny1) dz..... a eee 
Zo Zo 
oO oO 


Finally by substitution in Eq. (20) one obtains 


f Eoc (1 — e412) 
Emv = ered OG a ( 


0120 

Eq. (22) relates the plate efficiency Eyy to the operating 
conditions }, the point efficiency Eg, the mixing factor 
6 and the length of the tray z, 

A knowledge of the mixing factor ( is necessary, t 
enable the use of Eq. (22) to predict plate ica 

An experimental study of £6 as a function of the flow 
rates, the liquid properties, the weir height and some 
aspects of tray geometry was carried out in this work. 


NR 
Nm 


Experimental procedure 


(a) Mixing Factors. The apparatus consisted essen- 
tially of a single perforated plate 36 in. long and 3 in. 
wide, with 4 in. perforations spaced on a @ in. triangular 
pitch. The free area was 9.2% of the total bubbling area. 
The plate was common to two flow circuits as shown in 
Figure 2. 

The first circuit, carrying the liquid, consisted of the 
perforated plate PP, the constant head tank CHT, the 
supply tank ST, the flow meter V3, the centrifugal pump 
WP and the necessary connections. 


The second circuit, carrying the gas, consisted of the 
plate PP, the air blower B, the damping tank DT, the 
rotameters V1 and V2 and the necessary connections. 


To measure splashing the holes of the tray were 
blanked except for a small section 3 in. by 3 in. in the 
center of the tray. 

A Lucite box 3 in. by 3 in. by 18 in. was used to 
collect the liquid splashed (and it was able to be mounted 
upstream or downstream) as shown in Figure 3. 


The plate was operated at steady state conditions and 
the splashed liquid was collected for a period of 10 min- 
utes every one inch distance from the end of the small 
bubbling section. 
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Figure 2 — Schematic diagram of the apparatus used. 


List of symbols — 


PP = perforated plate 

LW = Lucite window 

CHT = constant head tank 

Ss = ammonia distribution coil 
ST = supply tank 

V1, V2 = valves to rotameters 

WP = centrifugal pump 

B = air blower 

M = eleciric motor 


DT = damping tank 

EC = exhaust chimney 

AC = ammonia cylinder 

£ = temperature measurement 
r = static pressure measurement 


The variables studied independently were gas and 
liquid flow rates, weir height, kinematic viscosity and 
tray geometry. The procedure was to keep all variables 
but one constant in turn. 

The viscosity of the liquid was changed by using 
sucrose solutions of different concentrations. In all other 
runs the system air-water was used. Viscosities and 
densities of the sucrose solutions were measured at 25°C 
using an Ostwald viscosimeter and a pycnometer re- 
spectively. 

(b) Ammonia desorption efficiencies. Plate effi- 
ciencies of ammonia desorption for the system air-water- 
ammonia were measured to test Eq. (22 

The procedure followed was to supply ammonia at a 
constant rate to the constant head tank and desorb am- 
monia on the plate. The whole bubbling area of the tray 
was used. 

Liquid samples were ae from the bottom of the 
plate, by means of | plastic tube samplers, spaced 
along the centre line of vt length of the tray, to study 
the concentration gradient. The liquid samples were 
analysed for ammonia by titration. 

E quilibrium values for the system were obtained from 
Perry “) and they were corrected for temperature 
changes by using the relation”); T log y = constant. 

Point efficiencies of the desorption of ammonia by air 
were obtained from data published in reference (7) for 
small bubble-caps of the University of Michigan column 
and they were assumed to be valid for a perforated plate 
at the same operating conditions. 

The vapor phase concentration gradient was calcu- 
lated from the liquid concentration gradient, using the 
relation y = mx Eog. 

The average value of y to be used in calculating Eyy 
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Figure 3 — Above—Sampling set up for splashing rate 
measurements. Below—Lucite box sample collector (plan 
and elevation). 


was obtained by graphical integration of the vapor con- 
centration curves. 

(c) Sodium chloride gradients. To study the 
effect of eddy diffusivity in the liquid layer, the system 
air-water was used, while a constant rate of saturated 
sodium chloride solution was injected at the outlet of the 
tray. Backward distribution of the salt was measured by 
Ww ithdrawing liquid samples in the same manner as in the 
ammonia runs. 

The samples were analysed by titration with silver 
nitrate. 


Sample calculations 


For the purpose of illustration some sample calcula- 
tions are presented in this section for the more important 
quantities calculated from the experimental data. 

(a) Mixing Factors. In Figure 4 the distribution of 
splashing with distance was shown for a few runs. 

The average fraction Q, or Qs which was splashed 
was obtained by graphical integration of the distribution 
curve. This area was represented by a block fixed at one 
corner at the edge of the bubbling section as indicated. 
The horizontal distance from the center of the bubbling 
section to the center of the block was taken as the 
average distance of splashing. 
For example for G = 3.52 ft./(sec.), L = 1980 lb./(hr.) 
and 1 in. weir, Q, = 0.045, Q, = 0.016, z, = 3.00 in., 
Z, == 2.75 mn. ang 
Qp-Z~ — Qy,Z), = 0.045 3.00 —-0.016 X 2.75 = 0.106 ir. 

(b) Gas residence time for the ammonia runs. 


mae . . . - a 
The gas residence time is defined as t; = , where a 
uc 


is the gas holdup on the tray, cu.ft./(sq.ft. of tray bub- 
bling area) and ug, the gas velocity, ft./(sec.) 

In Table 3, first row of figures, for example, 
Pressure drop liquid holdup = 1.48 in. 
Froth height = 4.00 in. 
Gas velocity = 3.45 ft./(sec.) 


48 
_ 4.00 — 1.48 _ 9961 sec. 


Hence t, 
G 3.45 x 12 


(c) Number of overall gas phase transfer units, 
Nog Nog was taken from data published in reference 8 
for the absorption and desorption of ammonia in the 
University of Michigan column. It was assumed that the 
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TABLE 1 


AIR-WATER-SYSTEM 


M1x1nG Factor VALUES AND AVERAGE SPLASHING RATES 


Avg. splashing 


Avg. splashing 


Avg. splashing 


Gas Liquid forward backward distance, in. Weir Froth Mixing 
Rate Rate height height Factor 
ft./(sec. ) Ib./(hr. ) in. in. B, in. 
lb./(hr.) | % of flow Ib./(hr. ) © of flow | Forward | Backward 
1 
Tray No. 1, Bubbling area 3” x 3”, 67 holes, 14” d. 
| 
3:52 | 1,980 99.00 4.5 31.65 1.6 3.00 2.75 1.0 3.00 0.106 
EY 1,860 96.75 iy 44.70 2.4 2.85 2:00 1.0 2.05 0.086 
3.52 1,680 67.20 4.0 30.25 1.8 2.80 2.60 1.0 auto 0.065 
3.52 1,470 73.50 5.0 51.50 Be 2.85 2.50 1.0 2.25 0.055 
3.De 1,080 69.25 6.4 67 .00 6.2 2.75 2.50 1.0 2.50 0.021 
3.52 900 72.00 8.0 72.00 8.0 2.75 2.50 1.0 2.50 0.020 
3.52 540 62.20 110 54.00 10.0 2.40 2.50 1.0 Bike 0.067 
1.48 1,980 13.85 0.7 1.90 0.1 2.50 2.50 1.0 2.00 0.015 
2.16 1,980 39.60 7.0 15.85 0.8 bite 2.50 1.0 2.50 0.035 
2.84 1,980 71.25 3.6 31.70 1.6 2.90 2.50 1.0 yy 0.064 
3.52 1,980 99 .00 5.0 31.65 1.6 2.85 2.60 1.0 3.00 0.101 
4.40 1,980 158.50 8.0 51.50 2.6 3.00 2.50 1.0 3.50 0.175 
2.56 900 17.10 1.9 15.30 A.@ AE 2.50 0.5 bake 0.010 
2.56 900 37.50 4.5 36.00 4.0 2.50 2.35 1.0 2.00 0.017 
2.56 900 81.00 9.0 63.00 7.0 Lite 2.85 2 3.00 0.027 
2.56 900 103.00 10.8 81.00 9.0 3.00 3.10 2.0 3.00 0.045 
Tray No. 2, Bubbling Area 3” x 3”, 18 holes, 1%” d. 
3.52 1,980 128.50 6.5 63.40 a8 3.25 2.35 1.0 3.00 0.136 
3.52 1, 605 112.50 7.0 57.80 3.6 3.00 2.35 1.0 2.50 0.125 
3.52 1,260 70.60 5.6 63.00 5.0 2.85 2.35 1.0 2.25 0.042 
3.52 900 50.40 5.6 54.00 6.0 3.00 2.50 1.0 2.00 0.018 
2.97 1,980 109.00 5.5 39.60 2.0 3.25 2.60 1.0 2:19 0.122 
2.16 1,980 79.25 4.0 29.70 Ls See 2.25 1.0 2.a0 0.096 
Dies 1,980 3.96 0.2 0.00 0.0 2.90 0.00 1.0 1.50 0.005 
Tray No. 3, Bubbling area 3” x 6”, 174 holes, 14” d. 

3.52 1,800 135.00 Ts 126.00 7.0 2.48 2.65 1.0 2.75 0.014 
3.52 1,080 91.75 8.5 91.75 8.5 2.75 2.65 1.0 2.75 0.009 
TABLE 2 
AIR-WATER-SUGAR SYSTEM 
MIXING Factor VALUES AND AVERAGE SPLASHING RATES (TRAY No. 1) 

Avg. splashing Avg. splashing Avg. splashing 
Gas Liquid Forward Backward distance, in. Weir Froth Liquid Mixing 
rate rate Height Height u/p factor B 
ft./(sec.) | Ib./(hr.) in. in. sq.ft./sec. in. 
Ib./(hr.)| © of flow Ib./(hr.)) © of flow Forward Backward at:25 c 
3.52 1,600 72.0 4.5 52.8 aa 3.25 3.50 1 3.00 0.13 x I0-* 0.036 
2.67 1,600 D1 .2 S.2 44.8 2.8 3.50 3.50 1 3.00 0.13 x 10-4 0.025 
2.16 1,600 41.6 2.6 35:2 Zon 3.50 3.50 1 2.40 0.13 x 10-* 0.014 
3.32 1,600 51.2 asa 41.6 2.6 3.50 3.50 1 2.50 0.36 x 10-4 0.021 
2.67 1,600 48.0 3.0 41.6 2.6 3.50 3.50 1 2.25 0.36 x 10-* 0.014 
2.16 1,600 5D.2 2.2 32.0 2.0 3.10 3.00 1 2.00 0.36 x 10-4 0.008 
3.52 1,600 76.8 1.8 75.0 4.4 3.350 3.25 1 2.350 0.26 x 107-4 0.025 
2.67 1,600 64.0 4.0 58.2 cit 3.50 ee 1 2.50 0.26 x 10-4 0.020 
2.16 1,600 38.4 2.4 38.4 2.4 3.29 3.00 1 2.00 0.26 x 1074 0.012 
gas phase was dispersed in the liquid by the small (14” d.) Therefore for the run mentioned previously 
bubble-caps of the Michigan column, to almost the same Nog = 18.4 X 0.061 = 1.12 
degree as by the sieve-tray used in this work for the (d) Overall gas phase point efficiency, E,,.. Use 


same operating conditions. In other words, the point 
efficiencies for the two trays were assumed the same. 
Linear variation of No, with t; was found to apply 
within the experimental error and the data 
related as: 


were cor- 


Noa 12t;, for ug 1 ft./sec. 
Nog = 15.9 tafor ug = 2 ft/sec. 
and Nog = 18.4t, for ug = 3 and 4 ft./sec. 
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of the Equation Ej, = 1 — e-Noc was made 

Hence, Egg = 1 — e712 = 0.674 or 67.4% 

(e) Vapor concentration gradient. The vapor 
concentration gradients presented in Figure 10 were cal- 
culated from the corresponding experimental liquid phase 
concentration gradients (Figures 8 and 9) by using the 
relation y = mxEo«, where m is the equilibrium constant. 

The average vapor composition yay, was calculated 
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Figure 4 — Distribution of splashing with distance. Liquid 
rate = 1,980 lb/hr. l-in. weir. 
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Figure 5 — Effect of gas and liquid velocity on mixing 
factor. 


Tray No. 1, 3” X 3”, 67 holes, 4%"d. 
_ Tray No. 2, 3” X 3”, 18 holes, 4d. 


by graphical integration of the vapor concentration pro- 
files. It was believed that this method of calculation 
would give better values than a material balance on the 
plate, since a slight error in the outlet liquid composition 
results in an error as large as 10% in the yayg value ob- 
tained from a material balance. 


Results and discussion 


The results are presented in Tables 1 to 4 and Figures 
4to ll. 
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Figure 6a — Effect of air velocity on the mixing factor 
Air-water system, liquid velocity = 0.257 ft/sec., 


l-in. weir, 6 oc W215 


Figure 6b — Effect of liquid velocity on the mixing 
factor. Air-water system, gas velocity = 3.52 ft/sec., 
l-in. weir, Soc mys-* 

Figure 6c — Effect of seal on the mixing factor. 
Air-water system, Ba h4.2 


Figure 6d — Effect of kinematic viscosity on the mixing 
factor. Air-water-sugar system, liquid velocity — 0.342 ft/ 
sec. l-in. weir. 

Parameter: gas velocity: 
+ ug = 3.69 ft./sec. 
@ u, = 2.67 ft./see. 
° ug = 2.16 ft./sec. 


—0.57 
ML 
pa (**) 


(a) Mixing factors. The variables thought to affect 
the mixing factor were studied independently and results 
are tabulated in Tables 1 and 2 and illustrated in Graphs 
5 and 6. 

A logarithmic plot of the function vs. each one of 
the variables gave satisfactory straight lines, with slopes 
as shown in Graph 6. 

The mixing factor data for tray No. 1 were cor- 
related with an average deviation + 14% by the empirical 
equation 


B = 0.688 x 1074 ug? ®.uL?- 385, h*2, (“:) haa . (23) 
L 


A plot of Eq. (23) is shown in Figure 7. It should be 
pointed out that this equation holds only for tray design 
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TABLE 3 
Data FOR De SORPTION OF AMMONIA BY AIR 





~ | - 








| Lewis 
cer. | te | Lewis | Env Emv 
Gas Liquid Weir | Froth Liquid Gas | Froth| mG | Eq. | Noc | Eoc Emv | Emv | Eq corrected 
rate rate height | height | Holdupt res. | density ~ (23) Ref. % Exp’! Eq. | (23) for end 
ft./(sec.)| Ib./(hr.) in. in. in. time | g/cm’ 4 in. (7) % (3) | & effects 
sec. % | % 
3.45 1,980 1 4.00 1.48 0.061 | 0.370 | 0.151 | 0.039 | 1.12 | 67.4 71.0 Tie) Face 
2.84 1,980 1 4.00 1.52 0.073 | 0.380 | 0.128 | 0.023 | 1.34 | 73.2 74.7 76.7 76.6 
aioe 1,980 1 3.50 1.38 0.076 | 0.392 | 0.104 | 0.019 | 1.20 | 70.0 73.0 12.3 12.6 
3.69 1,050 2 6.50 1.97 0.102 | 0.303 | 0.316 | 0.096 | 1.88 | 84.9 88.5 27.5 97.4 
3.45 1,050 2 6.00 Lad 0.102 | 0.295 | 0.284 | 0.084 | 1.88 | 84.9 96.3 95.7 95.6 - 
2.84 1,050 2 6.00 2.00 0.117 | 0.335 | 0.240 | 0.058 | 2.16 | 88.5 94.6 98.5 98.5 
3.69 1,050 3 7.00 Zan 0.102 | 0.360 | 0.322 | 0.310 | 1.87 | 84.6 96.0 97.2) O.i 
3.45 1,050 3 6.50 2.48 0.097 | 0.382 | 0.290 | 0.278 | 1.79 | 83.2 96.0 | 94.3 94.2 
2.84 1,050 3 6.50 2.48 0.118 | 0.382 | 0.264 | 0.184 | 2.17 | 88.6 | 102.5 | 100.0 | 99.8 
2.44 1,050 3 6.00 2.48 0.120 | 0.413 | 0.229 | 0.133 | 1.92 | 85.3 91.0 | 94.3 94.2 
36.9 560 4 7.50 2.64 0.110 | 0.351 | 0.537 | 0.218 | 2.02 | 86.8 | 98.7 | 110.0 | 110.0 99.5 
3.45 560 4 7:25 2.56 0.113 | 0.353 | 0.484 | 0.208 | 2.09 | 87.6 | 99.3 | 109.0 | 109.0 98.5 
2.84 560 4 7.00 2.60 0.129 | 0.371 | 0.408 | 0.129 | 2.38 | 90.7 99.8 | 110.0 | 110.0 99.3 
2.44 560 4 6.75 2.64 0.140 | 0.391 | 0.353 | 0.090 | 2.23 | 89.2 96.7 | 105.0 | 105.0 95.0 
tCalculated as the ito: of wet and dry pressure e drop. 
TABLE 4 
SUMMARY OF F CALCULATIONS FOR PRED DICTION OF F ME THYL- [soBuT TYL- Ket TONE EFFICIE NCIES Data OF ‘RE FERENCE “7 
Gas rate, ft./(sec.)..... 6.13 6.07 6.23 6.13 6.08 6.13 6.23 3.19 gi23 
Liquid rate, gal./(min.). ; 10.05 2% 05 10.05 10.05 38 38 38 38 38 
Weir height, in... .. jt a > 5 1 3 5 1 3 
Liquid residence time, sec......... be ite 7.61 10. 52 12.28 12.28 12.28 12.28 4.90 3.50 4.87 
mG 
i = \ bau Sikes |, eae 2.29 2.35 2.42 0.625 0.641 0.694 0.305 0.343 
Length of liquid path, in...... 15.6 15.6 15.6 186 | 156 15.6 15.6 15.6 15.6 
se x 10 sq.ft./(sec.).. : ..| 0.108 0.108 0.108 0.108 0.108 0.108 0.108 0.108 0. 
PL 
Mixing factor 8, in. (eq. 23). 0.042 1.935 | 121.0 118.0 0.424 | 16.7 102.0 0.061 2.91 
Eog, ref. (8), %.. 30.5 41.5 44 43.5 33 43 49 45 58 
Env, (complete mixing ) fi inca 30.5 41.5 44 43.5 33 43 49 45 58 
Emv, (Lewis no mixing) % 44 69 77 77 36.5 49 59 48.5 64 
Emv, (Exp’l ref. 8) % ‘ 47 64 60 65 | 34 | 45 56 46 50 
Emv, some mixing predic ted by Eq. (22) 
in this work, 9 Be is sees etre asi 44 65.7 91.5 51.3 36.5 | 48 53.2 48.2 63.8 


No. 1 having 4 
pitch. 
The range of the variables studied was: 


= 0.115 to 0.423 ft./ 


in. orifices on @ in. equilateral triangular 


= 1.48 to 4.4 ft./(sec.), ur (sec. ) 


h = 0.5 to 2.5 in., = 0.13 x 10‘ to 0.36 x 10~* sq.ft./ 


pL 


(sec. ) 


Qualitatively it could be stated that by decreasing the 
free area of the bubbling section splashing is — 
while no effect of the wall was observed for trays No. 
and No. 3 (Table 1). 

However, it was felt that more data are required 
before a quantitative conclusion could be drawn, regard- 
ing the effect of tray geometry. 


(b) Plate efficiencies. Plate efficiencies and liquid 
concentration gradients were measured for the system 
ammonia-air-water as a function of the gas rates, liquid 
rates and weir heights. These efficiencies checked within 
+ 2.3% with the values predicted by both the Lewis 
Equation and Eq. (22) proposed here, and are presented 
in Table 3. 


On that basis the conclusion was drawn that in these 
runs, the splashing rate was too small to have any effect 
on the Lewis efficiency. 
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For the 4-in. weir (Figure 9) the concentration grad- 
ients appeared to level off towards the end of the tray. 

This was thought to be due to end-effects of the high 
weir, causing some mixing of the liquid at the end of 
the tray. A plate length correction was applied to the 
Lewis Equation as suggested by Stone “*); for this par- 
ticular case the tray was taken as 0.25 foot shorter and 
the results were correlated fairly well. The dashed lines 
in Figure 9 represented the concentration gradients cor- 
rected for weir effects. 

In addition to the ammonia efficiencies studied, Eq. 
(17) was tested with data presented in reference (7) 
for the desorption of the methyl-isobutyl-ketone from 
water, (Table 4). Predicted efficiencies agreed within 
+7.0% to experimental, compared to + 10% for the 
Lewis prediction and this led to the conclusion, that 
mixing of the liquid by splashing may be the predomin- 
ant mechanism for a bubble-cap plate as well. 

However, splashing data are not available for bubble- 
cap plates to compare with the data reported here for 
perforated plates and permit a direct test of the postu- 
lation. 

The 6 values for bubble-caps were calculated by 
Eq. (23) and although in a few cases wide extrapolation 


was used, agreement with the experimental efficiencies 


was satisfactory. 
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ammonia system, 


(c) Salt concentration gradients. It was observed 
that in cases where backward splashing was significant 
(6-10% of the liquid rate for low liquid velocities) then 
an upstream salt gradient (Fi ‘igure 11) was appearing even 
20 in. from the injection point; with no splashing present 


backwards (0-1% of the liquid rate for high liquid 
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Figure 10 — Calculated vapor phase concentration 
gradients for ammonia system. 


velocities), no salt diffusion was found in the liquid layer 
beyond 2 in. from the injection point. 

It was concluded that the backward distribution of 
salt was caused mainly by splashing and that the eddy 
diffusivity in the liquid layer had no effect on back 
diffusion. The salt concentration gradients were pre- 
dicted fairly well by the relationship 


O 
a (24) 


Qrzr 
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Figure 11 — Concentration gradients for backward 
mixing of NaCl. Gas velocity 3.45 ft./sec., l-in. weir. 


which was derived from splashing considerations by the 
use of simplifying approximations and should therefore 
be considered as semi-theoretical. The constants k, and 
k, are defined by the boundary conditions. 


Conclusions 

With the apparatus used, the variables studied and 
within the range of the variables investigated, mixing on 
the perforated plate, as defined by the mixing factor 6, 
varied according to the empirical relation, 


0.57 
B = 0.688 x 10+ uoteneannes (HE) RPSL ENS (23) 
L 


Mixing factors measured ranged from 0 to 0.175 in., 
but theoretically may attain much higher values. Some 
relevant calculations are shown in Table 4. 

Mining on the plate was shown to be caused mainly 
by “jumping” of the liquid, while eddy diffusivity or 
mixing in the liquid layer at gas rates too low for measur- 
able splashing was shown to have no effect on backward 
diffusion (salt gradients). 


Eq. (22) and Eq. (23) are recommended to be used 
for the prediction of plate efficiencies, in the case where 
complete mixing in the vertical direction exists and 
splashing is significant. 

It was felt that a criterion for the use of Eq. (22) 
would be f values above unity, because only then Eq. (3) 
and (22) differ significantly. Plate efficiencies for the 
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desorption of ammonia from water, using air as desorb- 
ing medium, were measured ranging from 71% to 
102.5%. 

The mixing factors prevailing in these runs were too 
small to have any effect on the plate efficiency, therefore 
the plate was found to behave like the Lewis model. 


The test of Eq. (22) with data of reference (8) 
pointed to the conclusion that splashing of the liquid 
may be the mixing mechanism for a bubble-cap plate 
as well, since experimental efficiencies were predicted, 
using the mixing model proposed here, with an average 
deviation of +7% compared to + 10% of the Lewis 
equation. 

f 


It was felt that further information on splashing rates 
from different tray configurations would be necessary to 
quantitatively establish the effect of tray design and work 
along this line is under way. 


Nomenclature 
E i NN 6655586 vise hs CRs REED % or fraction 
G NE gcse au'gigth ss 4 wraee acne Sal Ib. moles/(hr.) 
h eI 0 5S ccs ngs seaveteaiord vss in. 
L AEDES oreo 6 oe Sink ham wetes Ib. moles/(hr. ) 
m = slope of equilibrium line 
N = number of transfer units 
Q = fraction of liquid rate splashing. Qr, downstream. 
Qs, upstream 
t = residence time on tray,............ sec. 
+ = temperature 
u = linear gas velocity based on total 
UMN UNE os 5 ise ore « oa eet ft./(sec. ) 
or linear liquid velocity based on 
PRCHARET CERO i523 secsis so 5 Seal aces ce 59 ft./(sec. ) 
Xx = mole fraction of more volatile component in the 
liquid phase 
y = mole fraction of more volatile component in the gas 
phase 
A = distance measured on the plate..... in. 
Zo = total length of tray bubbling area. . .in. 
Greek Letters 
a = gas holdup on the tray............ cu.ft./(sq. ft.) 
B oe EUROS SOI oo ina as sass 4 in. 
Y = activity coefficient. 
mG 
x = 
i 
iM WINES ss tre loe Ps ae id eae elo oo lb./(ft.) (sec.) 
p INE Koes 0 6 hk SID CRE Ib./(cu.ft.) 
L/p = kinematic Viscosity..........sce008 sq.ft./(sec. ) 
Subscripts 
avg. = average 
B = backward 
F = forward 
G = gas phase 
i = inlet 
L = liquid phase 
MV = Murphree vapor 
n,n +1 = platen, n + 1 
oO = overflow from tray 
OG = overall gas 
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High Purity Oxygen For 
Steel Making’ 


J. T. HUGILL? 


Recent developments in the manufacture of steel 
require the use of large quantities of high purity 
oxygen. To supply this demand new air separation 
plants have had to be designed and built which are 
capable of producing oxygen of large tonnage quan- 
tities at a minimum purity of 99.5%. This paper will 
describe in detail the 100 ton/day oxygen units built 
for Dominion Foundries and Steel Company, at 
Hamilton, two of which are now in operation. The 
process flowsheet and the reasons for adopting it is 
explained. Some of the problems encountered in the 
choice of certain specialized pieces of equipment are 
considered. Details of flow, pressure and tempera- 
tures are presented for normal operating conditions 
as well as operating control, methods, and measure- 
ment. The paper is illustrated with pictures of the 
equipment and the assembled plant. 


HE use of relatively pure oxygen in the manufacture 
7a steel has been tried in many ways over the past 
25 years. The first attempts were in Germany with a 
Bessemer convertor using an oxygen enriched blast. This 
showed considerable promise in increasing production 
but the advantages were offset by higher refractory con- 
sumption. On this continent the use of oxygen did not 
start until after World War Il. Many successful tests 
were run using oxygen accelerated combustion in oil or 
gas fired burners of the open hearth furnace and direct 
oxygen lancing of the open hearth charge to reduce 
carbon content rapidly. Some tests were done on enrich- 
ing blast furnace air with oxygen and one large oxygen 
plant was built to supply a ferro-manganese blast furnace 
with oxygen for continuous production. While many of 
the tests were successful in that the production was in- 
creased by a marked degree, the steel industry as a whole, 
until recently, has been slow in putting oxygen pro- 
cedures on a production basis. One process developed in 
Austria has however caught on faster than any others. 
This consists of top blowing with a jet of high purity 
oxygen into a charge of molten pig iron and producing 
a relatively high grade steel directly. This was the pro- 
cess decided on by Dominion Foundries and Steel Com- 
pany (Dofasco) after a very thorough investigation of 
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all possible methods for increasing steel production. The 
major advantage of the oxygen jet process of top blow- 
ing is that the basic capital cost is considerably lower. 
For the same tonnage output of a 250 ton/day open 
hearth furnace, the capital outlay, including the oxygen 
plant, is halved and the operating costs are somewhat 
lower “), From all reports the steel is of a superior quality. 


Oxygen plant 

The original plants using an oxygen jet process were 
based on 95% purity oxygen. It was determined that 
higher purity oxygen gave better steel due to the reduced 
nitrogen content. For example, increasing the oxygen 
purity from 99.2 to 99.5% reduced the nitrogen content 
in steel from 0.0035 to 0.0024 or a 30% reduction. This 
is particularly important in cold rolling light gauge sheet. 
Dofasco specified an oxygen purity of 99.5% or better. 

The next problems to be studied were oxygen pres- 
sure and its use in the steel plant. These two had to be 
considered together. The steel furnace operates on a 
batch basis for a total heat time of from 45 to 55 minutes 
with oxygen blowing taking from 20 to 25 minutes dur- 
ing each period. The rate of oxygen consumption is 
about 1800 scfm per ton of charge at an operating pres- 
sure of 125 psig. A low temperature oxygen plant, how- 
ever, operates on a continuous basis and while the average 
tonnage plant can operate between 70 and 110% of de- 
signed capacity, it takes up to one day to change the 
rate and, at the same time, keep the purity steady, 
especially if it must be maintained at 99.5 + 0.2%. The 
furnace consumption corresponds to a rate of 250 tons/ 
day for 20 to 25 minutes period in each hour but the 
total daily requirement would be only 80 to 85 tons. 


Based on this demand it was decided to design an 
oxygen plant to give a daily production of 100 tons of 
gaseous oxygen which would be compressed into a 
storage system at sufficiently high pressure and of suffi- 
cient volume so that the oxygen could be withdrawn at 
the furnace pressure for the peak flows. This was felt to 
be the best system for tying together a steady oxygen 
production rate with the intermittent high volume use. 


Oxygen compression and storage 

Compression of oxygen to high pressures of the order 
of 3,000 psig has been standard in the compressed gas 
industry for many years and provided certain design 
features are maintained, no serious troubles have been 
experienced. However, the volume involved has been 
quite small and only recently have the compressor manu- 
facturers been required to design machines capable of 
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compressing 1,500 scfm, or more, 


pressures in excess 
of 150 psig. The major problem in oxygen compression, 
using reciprocating compressors, is that absolutely no oil 
can be used. Water lubrication is standard practice in 
the compressed gas industry w hich, of course, entails a 
drier following the compressor. Carbon ring compressors 
have been developed to handle up to 2,500 scfm but 


water lubrication is necessary otherwise ring wear is 
excessive and consequently down time increases. Dry 
oxygen was required because of the corrosion problems 
resulting from the storage of oxygen in the steel storage 
vessels. If water lubrication was a necessity, a drier would 
have to be added to the circuit and it seemed rather use- 
less to take bone dry oxygen, add water and then remove 
it again, which removal would not be absolutely com- 
plete. Also the cost would be increased by the cost of 
the drier. 

The Sulzer Compressor Company had developed an 
oil free reciprocating compressor w ith a special labyrinth 
seal design and several smaller machines of this type were 
operating throughout Europe and giving very satisfac- 
tory service. They undertook to design and build a larger 
machine of the oil free type to handle 1,700 scfm to a 
pressure of 325 psig. Because continuous operation of 
up to one year was required, it was decided to install 
two full capacity machines. The storage sy stem had to 
be designed to handle peak flows at the maximum oxygen 
jet line pressure, allowing for control valves and line 
drop. This was specified at 175 psig and the volume re- 
quired above this pressure was 300,000 cu. ft. The stor- 
age system was therefore designed for a maximum pres- 
sure of 350 psig with a total capacity of 600,000 standard 
cu. ft. This consisted of 4 vessels each 10 feet in diameter 
and 70 feet high (Figure 1). 





Figure 1 — The first 100 ton/day oxygen plant at 
Dominion Foundries and Steel Company, Hamilton, Ont. 


Design of oxygen plant 
Several factors influenced the choice of cycle and the 

final design of the low temperature equipment: 

1. Low operating cost was essential. This meant low 

power consumption. The power consumption exclu- 

sive of oxygen compression was guaranteed not to 
exceed 400 kwh/ton of pure oxygen. 

Rapid delivery was specified by the customer. The 

plant had to be ready to go on stream in 13 months. 

This did not allow much time to consider new cycles 

and new equipment. 

3. The plant must be able to run for a minimum of 12 
months without stopping for defrosting. The equip- 
ment had to be very well proved in previous opera- 
tions. 

4. The oxygen purity must be 99.5%, or better. This 
automatically eliminated the standard regenerator 
type plant. 


Nm 


As a result of the above factors the basic cycle chosen 
was a modified regenerator type plant with a large 
booster cycle operating at a pressure of 450 psig. This 
allowed the use of regenerators to be restricted to the 
nitrcgen stream only. The expansion engine had to be 
a reciprocating type as at that time no turbo-expanders 
were developed to handle a pressure drop of 400 psi, at 
a relatively low specific volume. 


Description of cycle 


Figure 2 shows a flow sheet of the complete plant. 
The main air compression is performed with a Brown- 
Boveri Turbo blower. This is a 6-stage isotherm ty pe 
compressor with stage inter-cooling. The blower is 
driven through a speed reducing gear by a 2,000 H.P. 
synchronous motor and compresses the air to 75 psig. 
The air is cooled and scrubbed with water in a Raschig 
ring packed scrubbing tower. After cooling, the air 
stream is split and approximately 25% goes to a booster 
compressor where the pressure is further raised to 450 
psig. This machine is a two-stage reciprocating Ingersoll- 
Rand compressor. 

The 75% portion of the air at 75 psig flows to one 
of a pair of regenerators. These are packed with a large 
number of spirally wound aluminum pancakes. Details 
of the designs and construction of these elements have 
been thoroughly described in the literature (*). While the 
ete design of the regenerator is 25 years old it has yet 

» be equalled for this application. 


The air stream is cooled in one regenerator in contact 
with the cold packing. In the other regenerator, the 
separated nitrogen fraction is cooling the packing. The 
flows of air and nitrogen are reversed at a pre-deter- 
mined time, approximately every two or three minutes. 
The reversing valves are double port air operated, the 
air supply being controlled by solenoid valves which, in 
turn, are electrically operated from a timer. To avoid 
the large pressure changes i in reversal, an equalizing v valve 
opens during each cycle to equalize the pressure in the 
vessels. The whole reversing cycle lasts about 4 seconds. 
At the bottom of the regenerators, a system of check 
valves, especially designed to have a very low pressure 
drop, controls the flow of air and nitrogen in their re- 
spective circuits, thus the cold air from the bottom of the 
regenerator flows directly to the high pressure column. 
The portion of air which has been boosted in pressure 
to 450 psig is scrubbed with a 10% caustic solution to 
remove the CO,. The caustic tower is a standard design 
Raschig ring packed vessel with a large storage capacity 
at the bottom. A centrifugal pump circulates the caustic 
solution from the sump to the spray ring at the top of 
the tower. The volume of solution in the system is such 
that it can be circulated for eight hours before being 
discarded. The 10% solution being considered _ 
when it is 65% carbonated. A spare caustic pump 1 
provided together with a caustic make- -up tank and in- 
jection pump, the latter being used to pump fresh solu- 
tion into the system. Following the caustic tower and 
separator is a vessel filled with lump caustic which acts 
as a dryer to dry the air to approximately —60°F dew- 
point. The air then flows to one of a pair of switch ex- 
changers where it is cooled to about —90°F counter- 
current to the outgoing oxygen product and a small 
stream of waste nitrogen. The air flows in the shell and 
the products inside the tubes. The water contained 
the air below —60°F dewpoint is deposited in the ex- 
changer on cooling and when the pressure drop across 
the exchanger begins to build up above 10 psig the ex- 
changers are switched. One of these warm exchangers is 
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shown in Figure 3. The usual practice is to switch on a 
two week basis regardless of pressure drops. When an 
exchanger is taken out of service it is immediately de- 
rimed (defrosted ) ready for the next period of operation. 
The air is now further cooled in a liquefier following 
the warm exchanger. This is a shell and tube exchanger 
of similar design to the warm exchanger. Cooling is 
done also by countercurrent flow of oxygen product and 
waste nitrogen. Part way down the liquefier, the major 
portion of this air stream is withdrawn to the expansion 
engine. The take offs (each provided with a butterfly 
valve) are so designed that the temperature of the air 
can be closely regulated. This is important in getting the 
best efficiency out of the expansion engine. The inlet 
temperature of the expander must be such that when the 
air is reduced in pressure to 70 psig, and flowing directly 
into the high pressure column, the temperature of the 
air is between 3°F and 5°F above the liquefaction point. 
This gives the maximum enthalpy difference without 
liquid formation (*), Expanding into the liquid region not 
only reduces the enthalpy drop but creates a serious 
mechanical problem with the presence of liquid in the 
cylinder of the expander. This temperature level also 
gives the correct amount of reboil heat to the bottom 
of the high pressure column. Braking of the expander 
and recovery of the energy output is accomplished by 
driving a wound rotor motor at slightly above synchron- 
ous speed thus regenerating power and pumping this 
power back into the 550 volt line. The power recovery 
from the expander is more than sufficient to compensate 
for that used in the auxiliaries in the plant, caustic pumps, 
water pumps, etc. Speed control is achieved by a variable 
resistance on the motor windings and by varying the 
inlet valve cam settings. 

The balance of the high pressure air, now reduced to 
approximately 5% in volume, is further cooled and par- 
tially liquefied in the lower end of the liquefier and 
withdrawn at the bottom. It is circulated through a sub- 
cooler coil in the sump of the high pressure column and 
flashed as a subcooled liquid to the 7th tray. This coil 
immersed in the column bottoms provides further reboil 
for proper operation and at the same time insures that 
the air will be totally liquid when expanded from 450 
psig to 70 psig. Introducing a liquid stream at this point 
increases the reflux ratio of the lower part of the column 
and further improves the separation. 

The high pressure column is operated to give a bot- 
tom liquid product containing about 40% oxygen and 
an overhead of pure nitrogen. The condensed overhead 
is withdrawn, subcooled by the gaseous nitrogen product, 
and flashed to the top of the low pressure column to 
provide the necessary reflux. The bottoms or rich liquid 
is passed through one of a pair of silica gel filters, also 
subcooled and flashed as the feed to the mid point of 
the low pressure column at 5 psig. Using the low pressure 
cold nitrogen product gas to subcool the liquids reduces 
the flashing when these are expanded from the high 
pressure column pressure of 70 psig to the low pressure 
column at 5 psig. 

The purpose of the silica gel filters is to remove ace- 
tylene and some hydrocarbons from the rich liquid. The 
presence of small quantities of solid acetylene in liquid 
oxygen creates a very serious explosion hazard and over 
the years has caused many serious accidents. Activated 
silica gel has the property of adsorbing a limited amount 
of acetylene and thus prev enting or at least reducing its 
possible concentration in the oxygen. Acetylene is pre- 
sent in the air stream entering the plant and concentrates 
in the rich liquid (high pressure column bottoms). As 
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this is the feed stream to the low pressure column it is 
further concentrated in the liquid oxygen product. The 
solubility limit of acetylene in the liquid oxygen is of the 
order of 3-5 ppm, (5) above this concentration solid 
acetylene precipitates out and a potential explosive con- 
dition is set up. It is essential therefore that the filters 
reduce the acetylene content to a point well below the 
above figure. Details of filter design and allowable con- 
centrations have been thoroughly covered by Kerry(+ ® 
and others‘), Accepted practice is to use 100 Ibs. of 
silica gel per 1,000 cu. ft. of air being treated assuming 
a filter change every four to five days and a maximum 
allowable concentration of 1 ppm of acetylene in the 
incoming air. This means that on a 100 ton/day ~~ 
each filter will contain about 950 Ibs. of silica gel. A 
special pellatized form is used which reduces the loss by 
attrition. In areas where the acetylene content is likely 
to rise to a high figure it is necessary to introduce addi- 
tional removal measures. Obviously ‘the best place is on 
the incoming air stream. Location of the air inlet in 
relation to the sources of contamination is very im- 
portant and this must be studied for each installation. 
Recently a catalyst system has been designed to remove 
acetylene and some unsaturated hydrocarbons in the air 
feed. This equipment shows some promise and is the 
best possible answer to operating plants in areas where 
contamination can rise above normal safe limits. 


Main vaporizer or reboiler 


In an air plant the overhead condenser of the high 
pressure column is at the same time the reboiler of the 
low pressure column. Gaseous nitrogen from the top 
tray of the high pressure column at 70 psig and —287°F 
is condensed on one side while liquid oxygen at the 
bottom of the low pressure column is vaporized at 5 psig 
and —292°F. In the small standard type of air plant the 
main vaporizer forms the bottom of the low pressure 
column and the top of the high pressure column. The 
liquid oxy gen collects in this vessel surrounding the tubes 
of the vaporizer, inside of which the gaseous nitrogen 
from the top tray of the high pressure column condenses. 
With the rapid increase in size of air units in the past 
few years, vessels of this type become more complicated 
in their mechanical design, the principal problem being 
the tube sheet. To keep the already small temperature 
difference between oxygen and nitrogen, there is a limit 
to the permissible static head of liquid oxygen hence the 
larger the unit the larger the diameter of the vaporizer. 
The problem of designing a vaporizer with 12,000 half- 
inch tube on a diameter of over five feet becomes fairly 
great. A radical departure was made for larger scale air 
plants by adopting a falling film type of main reboiler. 
This vessel has a relatively smaller number of very long 
large diameter tubes down which the oxygen falls. 
Around these tubes the gaseous nitrogen condenses. This 
design however introduces another problem. In the stand- 
ard vaporizer, the low pressure column is mounted im- 
mediately above the high pressure column. In this way 
the oxygen flows by gravity to the main vaporizer. 
Placing the new design of main vaporizer between these 
two columns would increase the cold box height, adding 
complicated structural problems and most important of 
all increasing the losses by radiation due to increasing the 
over all surface of the cold box. It was decided to hang 
the main vaporizer tube nest in the top of the high 
pressure column (see Figure 2) thus giving a clear pas- 
sage up around the tubes for the nitrogen gas. By pump- 
ing the liquid oxygen from the bottom of the low pres- 
sure column to the top of the vaporizer it was possible 
to locate the LP column besides the HP column. The 
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Figure 3 — One of the completed warm exchangers in 
L’Air Liquide shop in Montreal ready for shipment to the 
job site. 


ee 





Figure 4 — The main control panel for the low 
temperature equipment. 


pumps are Lawrence centrifugal type with specially de- 
signed stuffing box, mechanical seal and long shaft so 
that the motor can be outside the insulation casing with 
the pump located sufficiently far inside to reduce heat 
losses. Because the fluid being pumped (liquid oxygen) 
is at a temperature of ~292°F at its boiling point some 
of the problems of packing and pump design required 
special attention. Two pumps are installed one in opera- 
tion and one spare. The spare pump is always kept 
cooled ready for immediate starting. Another very im- 
portant aspect of this reboiler design is the safety 
feature (*), The liquid oxygen collecting in the bottom 
of the colme is transferred continuously to the main 
vaporizer thus there is no possibility of a static bath of 
liquid where acetylene could precipitate out locally and 
accumulate. The rate of pumping is such that liquid 
oxygen is in continuous circulation. 

The oxygen product is withdrawn as a gas from the 
main vaporizer, passes through the tubes of the liquefier 
and warm exchanger where its cold is given up to cool 
the incoming air. ‘From the warm exchanger it flows to 
the suction side of the oxygen compressor. 

Reference has been made several times to the small 
stream of waste nitrogen used in the shell and tube ex- 
changer system together with the oxygen product. The 
purpose of this stream is two fold. The volume of oxygen 
produced is not quite sufficient to cool the high pressure 
air stream over a reasonable pressure variation. In the 


event that it was possible to operate this cycle at lower 





Figure 5 — No. 1 100 ton unit showing air compressor, 
regenerators and part of the cold box. 


pressure due to higher expander efficiency or lower than 
calcluated heat losses, there would be danger of a cross 
in one or other of the exchangers. With provision made 

handle this nitrogen stream and some minor modifi- 
cations to the column it would be possible to have this 
portion as relatively pure nitrogen, i.e. with an oxygen 
content of less than 1,000 ppm. 

Control of the oxygen suction pressure is maintained 
by an automatic pressure control valve — holding the 
pressure constant at 1.5 psig. This stabilizes the column 
Operation and assures a positive pressure to the com- 
pressor inlet. The compressors are equipped with tem- 
perature controllers on the discharge of each stage which 
act on the main switch to shut the machine down should 
the temperature rise above 460°F. This is considered the 
safe temperature limit for oxygen compressors. All pip- 
ing is copper and sized to keep the gas velocity at a 
low level. 

All the low temperature equipment is enclosed in a 
tight metal box structure commonly called a cold box. 
This consists of structural steel framing members on 
which are placed thin panels fastened w ith batten bars. 
The casing is fabricated to design dimensions by the 
Lindsay Canaan The box is completely packed with 
a special form of mineral wool to provide the necessary 
insulation. The depth of mineral wool between vessels 
and cold box casing dictates the size of the box, normally 
accepted practice is to have a minimum thickness of 
18 in. As the radiation heat losses are directly propor- 
tional to the square of the cold box surface the equip- 
ment is installed as closely as possible and hence the 
piping is somewhat complicated. Construction companies 
used to pipe racks and other standard methods of joining 
relating pieces of equipment would consider the piping 
in a cold box a “nightmare”. To further reduce the heat 
losses where the operating temperature is near or below 
air liquefaction temperature, some vessels have special 
foam glass insulation applied directly. The density of 
packing is very important in order to prevent heavy 
convection currents inside the box. 

All connecting piping is stainless steel and, except 
for flanged control valves, is all welded. Stainless steel 
meets the Code requirements for low temperature, it is 
easily obtained in all standard sizes and can be erected 
in the plant by the average qualified welder. All pressure 
vessels are designed according to the ASME Code re- 
quirements where they apply and are approved and 
stamped by the Inspectors of the Department of Labor, 
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or National Board before leaving the construction shop. 
This of course is standard practice in chemical industry 
as a whole but has only been adapted to low temperature 
equipment since 1947. 


Plant operation and control 


All the principal operating valves for the low tem- 
perature unit are air operated from the control panel. 
Automatic level recorder control is used on the rich 
liquid at the bottom of the high pressure column. With 
the big increase in the number of low temperature plants 
being built over the past few years it is now possible to 
obtain suitable valves and instruments designed for this 
service. Many other developments in the instrument field 
such as DP cells have been adapted to low temperature 
use with excellent results 7), 

Figure 4 shows a picture of the control panel for one 
of the 100 tons/day units. As will be noted most of the 
instruments are of the conventional type. A temperature 
recorder is used to record the important control points, 
namely the expander inlet and outlet and the regenerator 
mid point temperature. A 48 point temperature indicator 
enables the operator to check the whole process and to 
record on his data sheet the necessary control informa- 
tion. The timing devices for the regenerator reversing 
valves are supplied in duplicate with a gang switch so 
that if trouble occurs on one circuit the operator can im- 
mediately switch to the spare timer with only a minor 
upset in the plant operation. 

Measurement of the oxygen purity is done with a 
Beckman analyzer and recorded continuously. The waste 
nitrogen and reflux nitrogen are analyzed on a special 
design of Orsat by the operator. Analy sis of impurities, 
especially acetylene is a very important control feature 
and must be done regularly and accurately. Two systems 
are used, one by the operator on a shift basis which is 
qualitative in nature but indicates clearly and easily when 
the acetylene concentration is near the hazardous lev el; 
the other technique being somewhat longer is carried 
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out by the control laboratory once per day. Both analyses 
are done on the liquid oxygen from the LP column. In 
addition, regular checks are made on the rich liquid in 
each filter prior to reversal. 

The plant requires two operators per shift — one on 
the control panel and one on the compressors and pumps, 
The operator on the control panel is responsible for the 
overall running of the plant. In addition there is a fore- 
man and a maintenance man on day shift. 


With proper operation and a close control of the 
caustic scrubbing and regenerator reversing cycles, plants 
of this type should run continuously for at least one 
year. The first 100 ton/day unit has operated for periods 
up to 14 months without a shut down. Very little main- 
tenance is required on the cold box proper. The oxygen 
pumps require repacking about twice a year and replace- 
ment of their teflon bushings about once every 18 months, 


Since Dominion Foundries pioneered this process on > 


the North American continent there have been five other 
plants built by steel companies in the United States. 
Dofacso has added another furnace and a second 100 
ton/day oxygen plant. Both oxygen units are running at 
full capacity supplying the necessary oxygen for the 
continuous operation of two furnaces. 


The author wishes to thank Dominion Foundries and 
Steel Company Limited and L’Air Liquide for permission 
to publish this paper. 
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Some Experiments on Orifice Sprays 


WILLIAM E. RANZ 


New approaches to studies of spray systems were 
demonstrated by a series of experiments involving 
liquid jets spraying in liquids. Attention was focused 
on those conditions where the stresses causing break- 
up arose from the inertia of the surroundings and 
where the induced motion of the surroundings con- 
trolled the characteristics of the spray zone in front 
of the orifice. Experiments and analysis involved 
principle of balanced stresses, dispersion angle, 
motion of induced phase, development of spray 
zone, and drop size. ' 


Ree injected into a gas or into another liquid 
through a single hole orifice breaks my" into a spray 
and induces motion in its surroundings. Geometrically 
the system is the simplest of spraying devices. Phy sically 
the system is a series of stable dynamical states, all re- 
quiring complex fluid mechanical analy sis. 

Some of the experiments reported here are concerned 
with the development of a model of a spray from a single 
hole orifice. Instead of being injected into a gas, the spray 
liquid is injected into another liquid. In this way the 
time scale of the breakup process is so lengthened that 
modern high speed photographic techniques can follow 
and record the details of the process. Also in this way 
one can accent the important effects of two physical 
properties of the surroundings, density and viscosity, 
which are often forgotten in experiments where ambient 
gas conditions are relatively constant. Of itself, the model 
is important toward an understanding of emulsification 
processes (1), 

All of the experiments reported here are concerned 
with conditions for “atomization” or “fully developed 
spraying”, where the stresses causing breakup arise from 
the inertia of the surroundings and where the induced 
motion of the surroundings controls the characteristics 
of the spray zone in front of the orifice. 


Characteristic spraying stresses 


The best way to analyze a dynamical system is to 
consider characteristic stresses in the system. Balanced 


1Manuscript received March 7, 1958. 

2Pennsylvania State University, University Park, Pennsylvania. 

This article is based on a paper which was presented at the Joint 
A.L.Ch.E.-C.1L.C. Chemical Engineering Conference which was held in 
Montreal, April 20-23, 1958. 


stresses mark the limits of various stable dynamic states. 
Unbalanced stresess define the regions of acceleration and 
deceleration. Table 1 lists the characteristic stresses for 
an orifice nozzle spray: 


Table 1. Characteristic Stresses for an Orifice 
Nozzle Spray 


Force per Unit Area of 


Type of Stress: Interface: 
Ambient Phase: 
Impact or Inertial Stress peV"/28. 
Shear Stress 2uzV/D, 
Phase Being Dispersed: 
Impact or Inertial Stress —_p, V*,/2g. 
Shear Stress 2u,V,/D, 
Surface Normal Stress: 4o/D, 


where for the characteristic viscous stress the velocity 
gradient is assumed to be V,/(D,/2) and for the surface 
normal stress the pertinent radius is assumed to be D,/2. 
Because they are measurable, the characteristic length is 
taken to be D, and the characteristic velocity is taken to 
be V,. It must be remembered, however, that analysis of 
the various dynamical states requires that actual lengths 
and relative velocities be related to the characteristic 
values. 

It is logical to define what is popularly called “atom- 
ization” as that spray condition where the impact or 
inertial stress of the ambient phase greatly exceeds the 
surface normal stress of the head of the liquid column 
issuing from the orifice. In terms of the characteristic 
stresses of Table 2, atomization occurs when 


paV¥*.Do/Sdee SF Ao... +. on bianate detarer ater Om 


The onset of this type of spray is noticed at about 

peV*,D,/80g. ~ 1 as a lateral dispersion of the drops in 
he spray brought about by a waving or whipping motion 
of the liquid column (see Figure 1). 


When p,V*,D,/80g. exceeds unity, surface waves 
with wave lengths shorter than D, appear. These are 
generated by the inertia of the relative motion of the 
ambient phase and have a wave length, 1, which can also 
be estimated on the basis of balanced stresses (?) 


Pi Jaen dal 


where U is the actual relative velocity between phases. 
These waves grow rapidly into scale-like films and liga- 
ments which, being anchored on the main liquid column, 
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[ Mercury Sprays in Water 
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(A—Ref. (6); B—Ref. (7); C—Ref. (2). 


are drawn out to a rather thin size before the surface 
normal stress becomes dominant and collapses the liquid 
mass into drops. 

Balanced stress arguments can also be applied to the 
breakup region of the spray where the surface area is at 
a maximum and the relative inertial stress of the ambient 
phase is just balanced by the surface normal stress trying 
to form drops. Under these conditions 


pgU?D/80g. = 1... sie gohan 
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en F AE 


Ad 








35 120 


1 cm.) 


earbon tetrachloride sprays in water. 


where D is of the order of magnitude of the diameter 
of the resulting drop and U is the relative velocity at the 
point in question. To expand this concept into a theory 
for drop size and drop size distribution there remain the 
difficult problems of properly distributing the relative 
velocities and liquid fluxes throughout the breakup region. 

So far it has been assumed that the viscous shear 
stresses* are negligible. When either the viscosity be- 
comes very large or the length scale becomes very ‘small, 
for exampie, a small filament or drop size, causing the 
viscous stress to become of the same order of magnitude 
as the impact and surface normal stress, then correction 
must be made for the effect of viscosity. An estimate of 
a correction which accounts for the viscosity of the dis- 
persed phase can be obtained from related dynamical 
systems (see Figure 2). Corrections for the viscosity of 
the continuous phase will be of the same order of magni- 
tude for the same values of the appropriate stress ratios. 


Dispersion 


A spray drop at the moment of formation has a for- 
ward velocity and a lateral velocity. The forward velo- 
city of each drop is dependent on the original forward 
velocity of the mass of liquid from which it was formed. 
The lateral velocity of each drop is dependent on the 
original lateral siete of the mass of liquid from which 
it was formed. A distribution of lateral velocities results 
in random flight directions and a spray cone. 

Lateral velocities are caused by vorticity generated in 
the flow ahead of the orifice and in the flow at the point 
of breakup. Swirl devices, bends, screens, and naturally 
generated turbulence in the confined flow are responsible 
for vorticity created ahead of the orifice. Wave forma- 
tion accounts for vorticity generated near the point of 
breakup. 

Some insight into the lateral velocities generated at 
the point of breakup and the order of magnitude of the 
minimum dispersion angle may be obtained only con- 
sidering Taylor’s theory for wind generated waves ©). If 
the growth. factor of ‘infinitely small waves is propor- 
tional to the growth factor of finite waves and these, in 

*Viscous stretching and contracting normal stresses are also involved. 
The characteristic stress is written in the same way, but the velocity 


gradient is now taken in the direction of flow instead of at right angles 
to it. 
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turn, a°e proportional to the lateral velocity of resulting 
drops, also randomly distributed, one can write a pro- 
portiona lity for some average open angle. Thus, 
tan 6 ~ (da/dt)/U = ama!/U 

where # is the angle of dispersion from the centerline of 
flow and a’ is the amplitude at which the wave ceases to 
accelerate in its rise. Since the only linear dimension of 
Taylor’s system is Am, and Ay, is a logical breakup distance, 
one can predict that 

tan @ = function (amAm/U) 
For order of magnitude, the simplest equation is 

tan 6 © (1/A) (amd PCs ass aencte) 
where Gm and j,, are given by Taylor’s theory and where 
(1/A) is a constant of proportionality * be determined 
experimentally. A(tan 6/42) (pi/p,)'/* vs. Tay lor’s viscos- 
ity parameter, (pr/pg)o°g. Ju? :U*, as predicted by this 
order of magnitude equation, is plotted in Figure 3. 


Model system 

Consider now a model spray system where a liquid 
is sprayed into another liquid of density 10° times that of 
the gas. So long as the model system and the spray system 
have the same value of p,V*,D,/80g,., so long as the vis- 
cous stresses are negligible, then the liquid-liquid system 
should be dynamically similar in the forward direction. 
The advantages of the model are immediately apparent. 
For the same orifice and interfacial tension, for the same 
type of spray, the velocity in the model can be of the 
order of 30 times smaller. Presumably, the breakup pro- 
cess occurs in the slow motion of a time scale magnified 
30 times. 

Similarity in spray dispersion angle is probably not 
achieved unless the model and prototype operate at the 
same density ratio p,/p:. Because this ratio is so much 
larger in the model system, the spread of the spray is 
much more pronounced, and this aspect of spray forma- 
tion is exaggerated in the model. 

Because of the wide differences in physical properties 
between the model and an ordinary spray, the model is 
of value in checking the simplified theory developed 
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above. It also provides data from which more sophisti- 
cate theories can evolve. Because of the many variables, 
however, complete modelling, that is, equality between 
model and prototype of the necessary number of inde- 
pendent dimensionless groupings of variables, cannot be 
obtained. 


Experimental equipment 

The experimental equipment consisted of a 18x18x36- 
in. transparent plastic water tank, open at the top, and 
various nozzle systems for injecting immiscible liquids 
heavier than water. Spray liquid was injected through 
single hole orifices with orifice tube length to diameter 
ratio of 3 to 1. The. falling spray was captured in a 12-in. 
diameter glass funnel at the bottom of the tank, run into 
a collection cup, and drawn off through a tube in the 
bottom of the tank. Various sampling devices, such as a 
patternator, an impact tube, and a spray sampling cell 
could be placed in the spray. Photographs were taken 
through the plastic walls of the tank. The background 
was a frosted glass on the opposite tank wall behind 
which a constant light source or high-speed light was 
placed. 

Obtaining a representative drop sample was the most 
difficult type of experiment to perform. Figure 4 shows 
one of the successful trapping devices used. The spray 
was directed downward between rails along which a tray 
of sampling cells and cover could be moved rapidly. The 
cover and cell tray, being open at both ends, sliced out a 
portion of the spray volume and stopped at the center 
of the spray. The trapped droplets, whose total number 
could be varied by varying the distance between cover 
and tray, then settled into the sampling cells. Such a size 
distribution sample is a volume sample and may not rep- 
resent the drop size distribution being produced per unit 
time by the nozzle. 





) caine 
) Stop : 
) | 
Figure 4 — Trapping device for sampling emulsions. 
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Figure 5 — Distribution of liquid and momentum flux; m 
mercury spray in water. Figure 6 — Induced flow in a mercury in water spray as re 
Data: orifice diameter, 0.040 in.; total flow rate, 4.3 em.3/ shown by dye trace—arrows show start of free dye stream. m 
sec.; maximum impact pressure, 34 cm. water; lines show ts 
normal distribution curves with dispersion constants of 
a Sine Ae Induced flow patterns re 
Figure 6 shows the direction and nature of induced © 
Appearance of the spray water flow as revealed by a dye trace. Induced flow : 
Figure 1 shows high speed photographs of various ae at — pe: angles to a nye cme a ith- " 
stages of the development of a mercury spray and a in the spray a trace " aay =e " irection ; the dl 
carbon tetrachloride spray in water. Because there is spray axis and is somewhat dispersed in a turbulent 
only a small range of velocities between drip and hydro- fashion. A large scale wave motion is also shown to be . 
dynamic effects, smooth columnar breakup is of small ae in the ee This .— to re generated Dy a 
consequence in these systems. The appearance of atom- y — ant acting as a density aria and by a large . 
ization conditions, however, is very much like that of a velocity gradient in the radial direction. . 
liquid spray in air. Because of the low interfacial tension Drop sise be 
of the carbon tetrachloride-water sprays, hydrodynamic PI lhe scale tiene india ail f. 
effects appear before drip conditions are overcome. wm pig Ii hg ™ ae ; er eames * a1 
Tanasawa *) has suggested that the criterion for the ea” bie a a ; a4 =a oe 1 — so ee hn 33 
same spray “appearance” is equal values of (p,gV*,/ Fas ” Fone 7 oe snl an a eng agin the IS 
og.) (p/pg)"**, a conclusion he reached empirically. Con- outer cages of te ee arger droplets near the center. "I 
sideration of the atomization criterion and theory of The largest droplets are slightly larger than the orifice ” 
spray angle (see discussion of Figure 3) suggests an ex- a oo eee pee sc = oa taken : 
planation for Tanasawa’s “jet number”. However, be- with oz ade — evice Ww ae oa ae cells were , 
cause of the numerous quantities involved, complete aan. ~ er i egy ence fic; = ee — : 
dynamical similarity may be difficult to achieve. Appar- — a: 7 _ a die sb ee wad on ” 
ently, the model sprays are best suited for a study of ttl meee dt aad hes -_ . ie Streeton - Spray Grop 
fully-developed, “atomized” sprays and their component ee ee ee eee en eae ee 4 
dynamical processes. _ Figure 8 shows the size distribution of the spray of Is 
Figure 7 as a log-probability plot of size versus cumula- tl 
Spray and momentum flux tive volume per cent. V 
Figure 5 shows an example of mercury and total Figure 9 shows a size distribution and photomicro- 
momentum flux patterns across a mercury in water spray graph of drops in a carbon-tetrachloride in water spray 
at 6 in. from a 0.040-in. orifice where the mercury was taken with the sampler of Figure 4. In this case the size v 
injected at an average velocity of 530 cm/sec. The distribution for a volume of spray appears to be uniform " 
momentum flux at the orifice was 3.1x10* gm.mass cm./ throughout a large portion of the spray. a 
sec.?; the integrated force on the 6-in. plane was 3.0x10* : b 
dyne; as expected, momentum is conserved. Development of the spray zone in front of an 
To a large extent, both the liquid and momentum flux orifice nozzle 
are distributed across the spray axis in a normal prob- These experimental results suggest a spreading mech- 
ability curve, the tails of the curves being at slightly anism based on randomly distributed sideways velocities 
higher values of the distributed quantity than would be imparted to the drops by the breakup process and by 
expected from a normal distribution. Furthermore, the vorticity generated in the confined flow before the ori- 
dispersion constant of the momentum flux distribution is fice. The cumulative wakes of the droplets add up to a ' 
less than that of the liquid flux. considerable induced secondary flow, and the velocity of ? 
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the sp:ay is slowed in the forward direction as more and 
more ‘nomentum is transferred to more and more in- 
duced flow. 

If, for example, it was assumed only water flowed 
against the impact tube used to obtain Figure 5, the ob- 
served impact pressures would indicate a total flow of 
460 cm.*/sec. of water. If the 4.3 cm.*/sec. of mercury 
flowing had the same “average” speeds as the water when 
it passe d the measuring plane, at least 90% of the 460 
cm.’/sec. represents the maximum amount of water 
which could be flowing. Also, if only water flowed, the 
maximum velocity of the water when it passed the 
measuring plane would be 260 cm./sec., almost half the 
velocity of the jet issuing from the orifice. 

Consideration of the relative velocity between the 
drop and its immediate surroundings at the moment the 
liquid mass assumes an undistorted drop shape “*) requires 
a relative Reynolds’ number of the order of less than 
102 to 10%. The corresponding stopping distance, again 
with respect to a stationary surrounding medium, would 
be of the order of 30 to 100 times the quantity D,p,/18p,. 
While acceleration effects should be taken into account 
when the ratio p,/p, becomes much larger than 10-*, 
these statements are roughly applicable for all spraying 
systems. For a 1,000 micron mercury drop in water, the 
maximum relative velocity is 10 to 100 cm./sec., and this 
relative velocity must be achieved in a region where the 
momentum average velocity is still nearly equal to the 
injection velocity. 

In most sprays it would appear that the maximum 
relative velocity is less than 20% of the “average” velo- 
city of the w hole spray, reas there would be fur- 
ther breakup. Thus, as a rough approximation, it is 
assumed that at any point in a spray zone the momentum 
average velocity of the continuous phase is effectively 
the same as that of the dispersed phase. 

If the flux of the dispersed phase and the total 
momentum are normally distributed across the spray axis, 
and if the foregoing assumption about equal velocities is 
made, it necessarily follows 1) that the velocities are also 
normally distributed, 2) that the flux and velocity of the 
induced phase is distributed with the same dispersion 
factor as the flux and velocity of the dispersed phase, 
and 3) that the dispersion factor of the momentum flux 


is equal to V2 times the dispersion factor of the flux 
and velocity distributions, That this simplification is not 
entirely true is immediately apparent since, for the condi- 
tions, the spray liquid density in any cross section would 
have to be constant. We know this density goes to zero 
outside the spray and is probably lower at ‘the edges than 
in the center. 

These assumptions, however, permit a rational analysis 
of the development of the spray zone, and this analy sis 
is supported, in general, by the experimental evidence. If 
the spray liquid flux, spray liquid velocity, and induced 
velocity are distributed as 


w/Wim = Vs/Vim = Vg/Vem = exp (— 62r?).........(5) 


where subscript m denotes the maximum value on the 
axis, where vz = v, is assumed, and where 


Wim = p:*,V,d°/4 from a mass balance, a momentum 


balance can be written as 
fee) co 
pi(rD?2,/4)V2y = | wana + | oetamn Liane deta (6) 
oO oO 


Here the volume of the drops is considered negligible in 
calculating the momentum of the gas, and the induced 
air flow is assumed to enter at right angles to the spray 
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Drops collected at edge of spray. 


Figure 7 — Photomicrographs of mercury drops collected 
at various distances from the axis of a mercury in 
water spray. 

Data: orifice diameter, 0.040 in.; flow rate, 4.0 cm.°/see. 


Seale 3%” = 500 microns. 


axis. Substitution of Eq. (5) into Eq. (6) results in an 
expression for decreasing velocities in the expanding jet. 


Vi + 32 (pp/ps) (1/82D%,) — 1 

Ven/Vo = — 8 (pe pr) (1/8*D3,) ee 
where 1/6 = ztan@. z is the distance downstream from 
the orifice, and @ is the angle from the spray axis at which 
the flux has decreased to 1/e times its maximum value. 
The liquid density of the spray, p,,, is obtained as the 
ratio W,/V;. 


(7) 
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Figure 8 — Size distribution of a mercury in water spray. 
Data: orifice diameter, 0.040 in.; flow rate, 4.0 cm.?/sec.; 


pV", D,/8o8, 8; pV? Dmax/808, = 12.6 


Pts = Ws/Vt = Wem/Vem = ps2 2/4 (Vgm/Vo)......- (8) 
According to Eqs. 7 and 8, p,, is not a function of r and 
decreases downstream, at first being nearly inversely pro- 
portional to z, finally approaching an inverse square re- 
lationship. If spray photographs show side views of the 
spray, a constant density of drop images per unit area 
measured at various z-distances along the axis is evidence 
of inverse proportionality to the first power of the down- 
stream distance. Decreasing density in the z-direction 
indicates inverse proportionality to some higher power 
than unity. 


If continuity is preserved in the continuous phase and 
if no forward momentum is allowed in the surrounding 
medium, the induced flow enters the spray at right angles 
to the flow and at a considerable rate since the volumetric 

rate through a cross section is Vgm/d*®. The radial in- 
duced velocity of the continuous phase can be obtained 
from the differential equation of continuity 
10d 0 
+ ittig) > =~ Vg 


as (Gap = D5... 


where u, is the radial velocity in question and where the 
boundary conditions must satisfy overall continuity. 

The radial velocity outside the spray zone can be 
estimated to be that arising from a line sink on the spray 
axis. This velocity will have the value 


z= (1 ‘2rr) d (VgmT7, '6?)/dz 


2V o (z/r) tan?@ 


i; + 32 (p, Ips) (22 tan? 6 /D2,) 
forr > z tan 8 


Dispersion angle 

The equations developed above for dispersion and 
penetration are not particularly satisfactory since they 
contain an unknown quantity, tan @, which is also a 
function of the independent variables. However, 6 can be 
measured more easily than velocity, and once it is estab- 
lished, either experimentally or theoretically, these equa- 
tions become useful. 

When the viscosity parameter, (p:/pg) 0° g’./u?,U*, of 
Figure 3 exceeds a value of about 10%, as it does in the 
cases of liquid-liquid sprays studied here, A(tan 6/47) 
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Figure 9 — Size distribution of a carbon tetrachloride in 
water spray. 
Data: orifice diameter, 0.040 in.; flow rate, 7.2 cm.°/see.; 


0, V*,D,/8o8, = 220; peV*Dinax/808, = 66 


(p:/pg)'® is expected to have a value of V 3/6 where A 
is a constant to be determined experimentally. The data 
of Figure 5 give A ~ 20. Data from carbon tetrachloride 
in water sprays support this figure for the dispersion 
angle as defined, but do not support the theory with re- 
spect to the visible total spray angle. The tangent of the 
apparent total spray angle of the carbon tetrachloride in 
water spray was about two times that of the mercury in 
water spray, not three times as predicted by the theory. 
In the first case flux was lower in the outer reaches of 
the spray and in the second case higher than a normal 
distribution curve would predict. The edges of the car- 
bon tetrachloride in water spray were better defined 
than those of the mercury in water spray (see Figure 1). 
The carbon tetrachloride-water spray also appeared to be 
more uniform in density. In every case the dispersion 
angle appeared to remain constant with increasing 
velocity. 

As a further test of spray angles, a water in air spray 
from a 0.020 in. orifice was investigated. At a V, of 
290 ft./sec., tan @ at a distance of 11.4 in. from the orifies 
was found to be 0.007, giving an A value of approxi- 
mately 18 for the particular conditions. Since such a 
spray has a very narrow spray angle and a long breakup 
distance (electrical contact between the nozzle and an 
impact probe could be maintained up to 2.5 in.), such 
close agreement with liquid in liquid sprays seems 
fortuitous. 


In all these examples, the effect of vorticity in the 
confined flow has been assumed small. Although pre- 
cautions were taken to keep the flow laminar, it is not 
known what portion of the spray angle can be attributed 
to nozzle flow conditions. 


Agreement in the dispersion angles of relatively non- 
viscous liquid jets cannot be taken as much support for 
the theory as a general theory. In the case of liquid-air 
sprays where the viscosity parameter was less than 10°, 
Schweitzer’s data‘) give an A value of the order of 3. 
Furthermore, the spray angle increased with increasing 
velocity, in contradiction of the theory. Such poor agree- 
ment should not be surprising. The theory does not ac- 
count for lateral velocities generated in the flow and in 
the nozzle (an effect easily demonstrated), for breakup 
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length, for successive peeling off of liquid layers at 
smaller relative velocities for low impaction efficiencies 
at the sampling points, and for breakup conditions at 
velocities less than those which could be characterized as 
atomizing conditions. The theory in its present state is 
only a rational theory, but it is capable of refinement and 


study. 


Drop size distribution 

Figure 8 suggests a mechanistic explanation of the 
drop size distrtibution of a spray from a single hole 
orifice. The outer reaches of the spray are subjected 
to “atomizing” conditions, the size of drops being dis- 
tributed, perhaps normally, about the survival diameter 
for the relative velocity encountered. Successive layers 
peeling off under less and less relative velocity show 
larger and larger sizes until finally the core liquid re- 
maining, which is not subjected to high enough velocities 
to atomize, breaks up in a more leisurely fashion under 
the action of surface tension and relatively small fluid 
mechanical disturbances. 


Figure 9 shows the size distribution of a volume of 
a spray produced under much more intense atomizing 
conditions. If spray velocities have developed according 
to Eqs. 5 and 7, this size distribution for a volume of 
the spray is also the size distribution for the spray pro- 
duced by the nozzle i in a giv en time. Again there is noted 
a tendency toward a “maximum size” which in this case 
is much smaller than the orifice diameter. A Weber num- 
ber based on the maximum or median drop size (see 
Eq. 3) is not a constant, further evidence that fully 
atomized conditions are probably never obtained. 

Doubts about the accuracy of size sampling do not 
permit much trust in a quantitative interpretation of such 
data as that shown in Figures 8 and 9. These data do, 
however, accent the complexity of this aspect of spray- 
ing. As in the case of other spray parameters, it is doubt- 
ful that any meaningful correlation of size distribution 
can be obtained by dimensional analysis alone. 
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Table of Symbols 


A = constant 

a = amplitude of surface wave 

D = drop diameter. Dmax = maximum drop size 

D. = orifice diameter 

ge = dimensional conversion constant for use of mass and force 
units. g = 1 (gm.mass/dyne) (cm./sec.?) = 32.2 (Ib. 
force/Ib.mass) (ft./sec.?) 

r = radial distance from spray axis 

ug = radial velocity of continuous phase 

U_ = relative velocity between phases 

V. = volumetric average velocity based on orifice area 

Vg = velocity of continuous phase. vgm = maximum value 

vt = velocity of discontinuous phase. vin = maximum value 

w: = flux of discontinuous phase. Wim = maximum value 

z = axial distance from orifice 

a = growth factor of wave amplitude. am = maximum value 

6 = dispersion contant of normal distribution curve 

@ = dispersion angle at which the spray liquid flux has decreased 
to a value 1/e times the maximum value on the axis 

\ = wavelength of surface wave. Am = most probable value 

Mg = viscosity of continuous phase 

bs = viscosity of discontinuous phase 

Pp. = density of continuous phase 

pr = density of discontinuous phase 

pts = spray liquid density at a point in the spray 


= interfacial tension 
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The Relationship Between the Capacity 


and Efficiency of Dewaxing Filters 


R. M. BUTLER? 


The throughput, and the sharpness of separation 
of the oil from the wax, are factors which must be 
balanced against each other in the design and opera- 
tion of a lubricating oil dewaxing plant. Changes in 
this balance can be made by varying such quantities 
as dilution ratio, speed of filter rotation, quantity of 
wash applied and wax content of the feed. The re- 
lations between these quantities and their effect on 
the performance of a plant are discussed. The treat- 
ment is largely a theoretical one and involves an ex- 
tension of conventional filtration theory to take 
account of the quantity of liquid remaining in the 
filter cake and the rate at which wash can be applied 
to displace it. The predicted interrelationships be- 
tween the different variables have been shown by 
families of curves plotted on coordinates of dewaxed 
pil output (capacity) and dewaxed oil yield (effi- 
ciency). From these diagrams it is possible to pre- 
dict the effect of changes in the operating conditions 
on the yield and capacity for a plant. 


Ox of the steps in the manufacture of lubricating oil 
is that of dewaxing. In most cases this operation is 
carried out by diluting the waxy oil with a solvent, chill- 
ing the mixture until the wax has crystallized and filter- 
ing the resultant slurry in a rotary vacuum filter. The 
dewaxed oil is obtained by distilling the solvent from the 
filtrate and the wax product i is obtained by removing the 
solvent from the filter cake. 

The wax produced in such a process normally con- 
tains appreciable quantities of oil as a result of the large 
fraction of liquid which remains in the filter cake. This 
loss of oil in the wax is undesirable since it represents a 
loss in dewaxed oil yield; it also makes the further pro- 
cessing of the wax more difficult. The loss of oil in the 
filter cake can be reduced by the use of more solvent. 
This solvent can be applied either as dilution to the feed 
going to the filter or as wash on to the filter. However, 
for any particular set of operating conditions, there is a 
maximum rate at which wash can be applied before it 
runs off the filter into the pan. Variations in the other 
conditions affect this maximum wash rate thus, for 
example: 

1. Increasing the feed dilution ratio increases the 

wash rate by making the cake liquid more easily 


displaced. 


N 


Increasing the filter rotational speed increases the 
wash rate by decreasing the cake thickness. 


1Manuscript received March 29, 1958. 
2Research department, Imperial Oil Limited, Sarnia, Ont. 
This article is based on a paper which was prese »nted at the 40th Annual 


Conference of The Chemical Institute of Canada, 
1958. 
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3. The wash rate is increased by any decrease in the 

wax content of the feed. 

The solvent to feed dilution ratio also has another 
effect in that it changes the filter rate which is obtainable 
and may thus change the capacity of the plant. In general, 
it is found that the rate at which oil passes through the 
filter increases with dilution up to a maximum and then 
decreases. 

It is the purpose of this paper to show the relation- 
ships which may be expected between the following 
variables: 

(a) Yield of dewaxed oil. 

(b) Filter rate. 

(c) Feed Dilution Ratio. 

(d) Wax Content of Feed. 

(e) Speed of Filter Rotation. 

(f) Quantity of Wash Applied. 
This will be done by combining the equation for washing 
previously developed () with standard filtration theory. 
The resulting theoretical relationships will show how the 
various independent operating filter variables may be 
expected to affect the separation efficiency and capacity 
of a dewaxing filter. 


Theory 


A. Simple filtration on a filter leaf—no washing, no 
drying 
This is the simplest case of filtration. The equations 
which will be developed are extended to the case of 
rotary filters in later sections. 
Let A’ be the filter area 
V, be the volume of oil in the filtrate after time 6 


Xx be the volume fraction of oil in the liquid part 
of the filter feed 


w be the ratio of wax to oil in the feed 

m be the ratio of liquids to solids in the cake 

u be the viscosity of the filtrate 

V be the volume of oil in the filter cake after 


time 6 
It should be noted that the rate at which oil enters a 
filter cake is not the same as the rate at which it leaves 
since some oil remains behind in the cake liquids. This 
fact has been neglected by previous writers on the subject 
of filtration in dewaxing plants“), 
ae will be proportional to 


(1) the area of filter A’ 


(2) the fraction of oil in the filtrate, x 
and from Poiseuille’s equation inv ersely propor- 
tional to 


(3) the viscosity of the filtrate, u 
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(4) the thickness of the wax cake. This in turn is pro- 
portional to the quantity of wax per unit area, 


(Vi + Vor 
A’ 
These may be combined to give, 
Ve: x A’)? 
Be eM  cscsssxens ce 
dO6 = uw (Vi + Vow 
The total volume of liquid remaining in the cake is 
’ ** i 
Vi" and also m (¥, + Vader" 
z= 
Hence, 
, Vo 
Vi ee Te ee ah a a (2) 
1 
_— 
Xmw 
Eq. (1) and Eq. (2) may be combined to give, 
dV. k? X (A’)? 
= =. ats UO. MWS caccnteciusnd (3) 
d@ 2 m Vow 
<2 . . . . 
where k is a proportionality factor which depends upon 


the size and type of wax crystals which are present. 
Equation (3) may be integrated to give, 


2 oe ‘on cA) 
A - - 


This equation describes the variation of the quantity 
of filtrate oil with dilution, wax content of feed, and filter- 
ing time. In examining the equation, it should be remem- 
bered that the viscosity, «, is a function of the concentra- 
tion of the oil, x. 

At very low dilutions, xmw_ is greater than 1, and 
Eq. (4) indicates that the volume of filtrate is an imaginary 

number. The physical significance of this is that with these 
very small dilutions there is not sufficient liquid present to 
saturate the cake and the equation is not applicable. The 
rate is zero. At values of xmw less than one, the minimum 
amount of liquid is exceeded and filtration is possible. As 


1-xmw 


the dilution is increased, increases faster than x 


decreases and the quantity of filtrate oil increases. This 
continues up to a point where the quantity of filtrate oil 
ina giv en time reaches a maximum, beyond this point the 
quantity of filtrate oil decreases and approaches zero at 
infinite dilution. 

The variation of viscosity with concentration is most 
readily represented by the followi ing equation which is due 
to Reeves(), 


p= BO — xpi ec a ee.. ae » CS) 
where § and y are constants which are diepcieinn upon the 
oil, solvent, ond filtering temperature. The following ex- 
perimental values have been obtained for the constants 
in this equation for a phenol treated S.A.E. 30 grade 
lubricating oil dissolved in a ketone type dewaxing 
solvent at 23°F. 
y =—3.24 
B.= 0.656 centipoise 


x 
A curve of 
u 


The yield of oil based on wax free feed is given by 
ae Ve 
Vi +Ve 


versus x is given in Figure 1. 


ae Dok s oo cee 


° This and the following equations are derived for a constant vacuum. 
°° Because V; is the volume of all in the cake and the ratio of total 
liquid to dissolved oil is 1-x. 

Because the total volume of solids (wax) is (Vi 
liquids to solids ratio is m. 


= + Vo)w and the 
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This equation shows that yield increases with dilution. 


B. Simple filtration on a rotary filter—no wash- 
ing, no drying 

The above equations may also be applied to a rotary 
filter in which the cake is neither washed nor dried. 


Let a be the active filter area at any instant 


A be the total filter area 
n be the filter speed in r.p.m. 


The following substitutions may be made in Eq. (4) 


A’=a 
Rites a l 
An 


Let V be the volume of filtrate oil per minute. Then, 


V Ve P 
A = - n. ras e2 3 Ce) 

q. (4) becomes, 
kg agi. — xmw ; (8) 


Eq. (6) is applicable unchanged, 
Y =1 — xmw.... (6) 


A series of calculations for the case of dewaxing an 
SAE 30 grade oil has been made and the results are 
shown graphically in Figure 2. 

This graph shows the effect of wax content of the 
feed and of dilution upon the filter rate (barrels of de- 
waxed oil per ft.2 per hour) and upon the yield (based 
upon wax free oil). The curves are drawn for a liquid- 
solids ratio of 10, a value which is of the same order as 
that found for filter cakes before they are dried. The 
value for the filter rate constant k has also been chosen 
to be of the same order as that found in practice. The 
curves are for simple filtration with no washing or drying 
and have been drawn to serve as an introduction to later 
curves which will be drawn for washed filters. There 
are three principles which should be noted. 

(1) For a given feed (i.e., given wax content), the yield 
of oil increases as the dilution is increased. The rate 
of production of dewaxed oil reaches a maximum 
value with dilution ratios in the range of 4 to 6. The 
filter rate maxima are relatively flat, however. 

(2) For a given dilution ratio, increasing the wax content 
of the feed decreases both the yield and the filter 
rate. The reduction in filter rate is due mainly to the 
increasing thickness of the wax cake and partly to 
the lower yield. The reduction in yield is the result 
of the larger volume of liquid w hich remains in the 
wax cake. 

(3) The dilution ratio which gives the maximum filter 
rate increases with the wax content of the feed. 


C. Washed filter leaf 

In earlier work) is was shown that the degree of 
washing obtained could be calculated from the following 
relationship: 

Volume of Oil in Washed Cake __ 1 Volume of Cake Liquids 
Volume of Oil in Unwashed Cake 4 Volume of Wash 
provided that the volume of wash was greater than half 
of the volume of cake liquids. In applying this equation 
to actual problems it is necessary to be able to predict the 
quantity of wash which may be applied in a given time. 


Eqs. (3) and (4) show that the instantaneous rate of 
filtration decreases during the filtration cycle because of 
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Figure 1 — Variation of x/y with oil concentration for Figure 2 — Theoretical effect of dilution and wax content 
SAE 30 oil dissolved in ketone solvent at 23°F. on yield and filter rate for simple filtration, n = 0.5 r.p.m.; 0 
a—0.4A; m= 10; k = 0.0041 B.cp %/ft?min %; Feed- r 
stock SAE : ‘aXy ist. saci Ss wax c s. 
the increasing resistance to flow which is exerted by the stock SAE 30 — dist. of eee Se Sea 
No wash, no drying. \ 
filter cake. The slowest filter rate is obtained at the end : 
of the filter cycle. The rate at which wash goes through oe , é ; 
Dividing Eq. (11) and Eq. (12) 
the cake will start off equal to this final filter rate and e Eq: t85) 1. ( 
wash is less than that of the filtrate. In order to simplify Vo® 2x Vo OF 
the calculation, the assumption will be made that the or 
washing rate remains constant at a value equal to the final 1V. Ow 
filter rate. This assumption will not be too far in error = 5 ee vee (13) ' 


for the following reasons: 
The viscosity of the filtrate is not very much greater _ This gives us the volume of wash to substitute in 
than that of the pure solvent (30-40% at the ‘most). Eq. c¥)- — 

(2) In practical applications the quantity of wash which The volume of cake liquids is given by, 


is applied is often comparatively small and the vis- Vv ' _ m(V, + V,)w. 
cosity does not decrease to that of the pure solvent. , may be eliminated from this expression by means of 


Ex riving 
(3) It is difficult to spray wash on to a cs cake in such |: (2), 8 8; 


a manner that the surface remains uniformly flooded. Velotae<é cite Hania Vo (14) 
; “ ; s = Ps uae eee 
As a result of this, it is probable that the washing 1 


rate is somewhat less than the maximum theoretical — 
one, 


Substituting Eq. (13) and Eq. (14) into Eq. (7) gives 


If it is assumed that the rate of washing is equal to 
5 the result, 


the final rate of filtration, then the following expression 


may be derived from Eq. (3). Volume of Oil in Washed Cake 7 6/8, (15) 


dQ k? 1 (A1)2 Volume of Oil in Unwashed Cake : ( e ') 


=  } (10) 


dé 2 p Vow xmw 


Hence the fraction of oil in feed remaining in cake 


If this is constant then, after washing is 

<1 (A')? (1 — xmw) 
2u Vow mxw 9 
: , p 98 w 
where @w is the time allowed for washing - = 


will then increase somewhat since the viscosity of the Q 11 4 


and Q is the volume of wash which can be applied. o > 6 1 
But, ~ \mxw 
? (1 — xmw) »)\2 ; 
Vo? = k? (Al)? ~ rere. aaa eee 6 porter CURA alias (16) 
be W 2(1 — mxw) 6, 
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where Op is the volume of oil in the washed cake and 
O is the volume of oil in the feed. The yield of oil may 
be calculated from this equation, 


D. Fully washed rotary filter 

Primary dewaxing filters are normally operated with 
as much wash as possible applied to the upper filter sur- 
face. The equations dev eloped previously can be readily 
adapted to this important case. 

The total rate at which oil enters a simple rotary 
filter is given by dividing the rate at which filtrate oil 
leaves (Eq. (8)) by the yield for a simple filter (Eq. 
(6)). This result is, 


Vi 
= gear, 1 (17) 
\ A Bw w(l — xmw) 


where V; represents this total rate. This value will be 
the same irrespective of whether the filter is washed or 
not. 
For a washed filter, the fraction of the oil in the feed 
which passes through the cloth is given by, 
a 


(mxw )? 
) 


from Eq. (16) (18) 
2(1 — mxw) 
where b is the washed filter area at any instant 


Hence the total filtrate rate, V, for a washed filter is 
obtained by multiplying Eq. (17) by Eq. (18) with the 


result, 
ae 
w me — mxw) 


20 — mxw) — (mxw)? b; ree wate tenc Gael 


al 


2(1 — mxw) 
It is also of interest to be able to calculate the cake 
thickness. This may be determined from the relation 


Vi + 1) 
t= = See ee ee .(20) 
A n 
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Figure 3 — The effect of dilution and wax content on yield 
and filter rate for a washed rotary filter. n = 0.5 r.p.m.; 


a= b = 0.4A; m = 10; K = 0.0041 B.rep 4/ft?'min\%; 


Feed SAE 30 distillates of various wax contents. 
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SAE 30 distillate. 


where t is the cake thickness. Eq. (20) is dimension- 


less. If ‘ is expressed in B/ft.* hr. and n is in r.p.m., the 
following dimensional equation may be obtained 
Viw(m +1) _. : 
= 1.12 A (Dimensional)........ (20a) 
n 


where t is in inches. 


Diagrams showing performance of washed 
rotary filters 


A. Effect of wax content of feed and of slurry dilu- 
tion. 

Figure 3 is similar to Figure 2 in that it shows the 

effect of dilution and feed wax content upon the filtrate 
rate and yield for a ty pical dewaxing operation. The 
conditions are the same in each case except for the fact 
that the filter in Figure 3 is washed. Several points should 
be noted: 

(1) For a given feed increasing the dilution has the fol- 

a effects: 

(a) The yield increases partly because the liquid re- 
maining in the cake after filtration contains a 
lower concentration of oil, and partly because 
more wash can be applied to displace it; 

(b) The quantity of wash which may be applied 
increases as a result of the decreased viscosity of 
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the filtrate, and also because of the decreased (5) The amount of dilution also increases in proportion 
thickness of the cake; to the feed rate. 
(c) The output of filtrate oil rises to a very flat 


; Another case which is of practical importance is the 
maximum and then decreases. 


effect of changing the filter speed while maintaining the 


(2) For a given total solvent circulation, an increase in solvent circulation rate at a fixed level. An increase jn 
the wax content of the feed results in a decrease in filter speed under these conditions has the following 
the filtrate rate, and a decrease in the yield of oil results: 
on a wax-free basis. This must also be accompanied (1) A considerable decrease in the yield of oil, i.e., more 
by a slight increase in the dilution ratio and a re- oil is left in the wax; 
duction in the quantity of wash which may be ap- 5 ; i 

: (2) A relatively small increase in the output of dewaxed 
plied. There is more oil left in the cake because oll: ; 
there is more cake and the amount of wash which : ; ; 
can be applied to this cake is even less than it was (3) A small decrease in the cake thickness; 
previously. (4) A decrease in the dilution ratio; 
A comparison of Figure 3 with Figure 2 shows the (5) A very slight decrease in the quantity of wash which 
considerable improvement in both yield and output can be applied. 

which can be obtained by washing. Thus, although it is possible to increase the output 

from a dewaxing filter by increasing the speed of rota- 

B. Effect of filter speed. tion, this output can only be obtained at the same yield 

In Figure 4, the effect of varying the dilution and if extra solvent circulation capacity is available. An in- 
filter speed is shown for a given feedstock. crease in filter speed at constant solvent circulation re- 

An increase in filter speed with a constant slurry dilu- _ - : ray an in output and a considerable 
tion ratio has the following effect: decrease in the yield. 

(1) The yield is unchanged, Diagrams similar to those shown in Figures 2, 3 and 

(2) The output of dewaxed oil is increased in proportion + can be drawn to show the effect of other variables such 
to the square root of the Gleer speed; as the ratio of w ashing to filtering area, the composition 
a s j : ' of the solvent, and the viscosity of the feed stock. 

(3) The cake thickness decreases in proportion to the 
increase of the square root of the filter speed; 

(4) The amount of wash required increases proportion- References 
‘ , > > ate: (1) Butler, R. M. and Tiedje, J. L., Can. J. Technol. 35:455-467, 1957. 
ally to the feed rate; (2) Reeves, E. J. Ind. Eng .Chem, 39:203, 1947. 

SUMMARY OF EQUATIONS 
Simple Filtration Filter Leaf 
Vo a V “peegeerne 
Oil through Filter _=k - 1 — xmw 6 ih sone 
A e A 
Yield Y =1-—-— mxw Y aad vee 
F F ees 
Waxy Oil to Filter ~=k(1+w) = 6 ; =k(1+w), /@ - 1 
A uMw(l — xmw) A A Nuw(1 — xmw) 
Fully Washed Filter 
. | | Os . a) 
Vi |x @ 2(1 — mxw) — (mxw)? Vi an x 2(1 — mxw) — (mxw)? =| 
Oil through Filter \ =k \ 1 ) Ow : \ =k be 
B Ww — xmw) . _ 
B ” { 2(1 — mxw) | pw ( — 2(1 — mxw) | 
: (mxw)? Or (mxw)? sa 
Yield Y -1- ii ie oh 
2(1 — mxw) 6y 2(1 - mxw) b 
‘ V i V 
Total Oil into Filter : =k |x eA an 
A V uw (1 — xmw) A uw ( — xmw) 
f ») k Ow | ee kb |; a gree 
Wash on to Filter Q == | Oy \- 1 — xmw Q = > fan [x1 xmw 
A 2x Or WV A 2x a A Nu Ww 
V Vi w (m + 1) 
Cake Thickness t =—w (m + 1) t =_—.- nis 
A A n 
ri s em. 1 
Dilution ratio D = ; or 3 es 
x (i ~ 1+(w4+1)D 


: Vi w 1 ee 
Dilution Solvent BS vec Z D ah aes 


A a. OF A xY 









¢% Oil in Wax : x 100% All of the above — oe 
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